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N ATTEMPT TO CALCULATE THE REACTION “N(d, p)5N* BY THE 
BORN FIRST APPROXIMATION 


By Ryszarp KOLODZIEJSKI 
Institute of Theoretical Physics, Warsaw 


(Received July 29, 1959; revised manuscript received October 26, 1960) 


The stripping differential cross sections are calculated for the case when the 25 level of 
15N is formed. Thomas' variational method is applied and the angular distribution of the emerg- 
ing protons is calculated by the Born first approximation for four values of the range of the 
Yukawa potential. The proton-target nucleus forces are found to be at least just as important 
as the neutron-target nucleus forces. The angular distributions as obtained for the emerging 
protons and the absolute values of the differential cross sections are entirely different from 
the experimental values. The calculations give some hint that in stripping reactions the capture 
process takes place only on the edge of the nucleus. 


The differential cross section for the reaction N(d, p)15N* was measured by War- 
burton and Ме Gruer (1956), who used 14.8 MeV deuterons and measured the stripping 
differential cross section also for the case when the 25 level (E —8.32 MeV) of the 5N nucleus 
is formed by neutron capture. It is interesting to compare their experimental results with 
the results of calculations carried out by the Born first approximation. In the present calcu- 
lations, both proton-neutron interaction (potential Vpn) and proton-target nucleus interac- 
tion (potential Vpr) are taken into account. Thus, in our calculations we retain the two terms 
<f\Vprli> and <f|Vpyli> in Owen and Fulton’s exact formula (1957) for the stripping 


differential cross section: 


"do а О) > (Ma M s 
dQ 2 QJ441) Сл): (2М+ Mr)? va 


+ Cf [Venlo i» + CfIAGH DP (1) 


KI Verl i> + 


Eq. (1) is identical with the Owen and Fulton (1957) formula for a target-nucleus of finite 
mass Му averaged over the initial spin orientations and summed over the final spin orienta- 
tions. The symbols J, and J therein denote the ground state spin of MN and the 25 level 
(E — 8.32 MeV) spin of 15N, respectively. The velocities v, and од are the relative velocities 
(as referred to the centre of mass of the final or initial nucleus) of the proton and deu- 


teron. 
(215) 
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The operators ША апд Н; in the third term of eq. (1) are given Бу 
Hi = Vyrt Vpr (2 


і a - PN + J PT (3: 


The symbol G in the third term in eq. (1) denotes the Green function for the full Hamilto: 
nian of our problem. The matrix elements of eq. (1) are taken between the initial and fi 
state wave functions. These are of the form 


y; = (2л) "^O (rz — ry)C(E) exp [i$ kan + Tp) — т). (4 
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yy = (22) ^^ (їм — rr) (E) exp и» ПАД, - M ikpfN Mi ла З (5 


. = > . . ET H us ^ © E сев 
The function Ф (ту — rp) in (4) is the wave function of the deuteron. | | д 
The functions С(Е) and ny(ry — 77) (5) are the ground state wave function of MN anc 
the 25-level wave function of 15N, respectively. The wave vectors ky and k, are the deuteror 
and emerging proton wave vectors in the centre of mass system. The a of these 


vectors are: 
pre ANT 2Мава \ | EY? Il 
4 Ма+ Мт h? са 


2Мт %„. 14.8MeV |^ 
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E | SEN ам 
kp = | h2 | E 216 ^ 10 cm “2.296 MeV А (7) 
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We assume 
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where М is the nucleon mass and M, — that of the deuteron. The energy E, of the emergin 


proton in the centre of mass system is given by 


MrEQO __ | .2M - Mr 
Mat Mz N= Mp, Be 8 


where e, = 2.226 MeV, and ey = 2.518 MeV. 

Eq and ey denote the deuteron binding energy and the binding energy of the captured neutron.: 
In order to calculate En We assume the value 10.834 MeV for the last neutron binding energy: 
in the ground state of 15N (Ajzenberg and Lauritsen 1952), 


For the relative velocities ГАСУ гапо, we have the formula 


Фр. (2M + Mr) Mr Makp "T 2M rk, 9) 
Vy (Ма + My) (M3 My) Mk, ^ (M+ Mr) ka (9) 
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‘hus, the formula for the Born stripping differential cross section is finally obtained in the 
orm 
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n order to evaluate the matrix elements <f|Vprli> and <f|Vpnli >, we assume, for the 
leuteron, the Hulthén function: 


(лн cl e Pr 


Фу» — in) = Na“ т , (11) 


стоп &G —2.316 x 10!? cm-!, and B=12.68 x 101? cm-1. 
The normalization constant N, is given by 


Nes E 2l (12) 


[The Yukawa potential / ру associated with this function has the depth P, and the range ro 


"ven by 
__ ET 4C 
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Therefore, the potential Vp, сап be written as follows 
ekalrp-"n| 


[Fp — ?м| 


(15) 


Vpn = ВРМ 
Ihe constant gp, which, similarly to the constant g of the Vpr potential, appears in our 
formula (10), can be expressed thus: 
K2 
SPN — па ? 


wherein 
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It is convenient to assume the Wilson function for the bound neutron. We have thus 


NN T 2 = Nye? NT, (16) 


with 
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Assuming the potential Vyr to be of the form 


і e "IrN-rTI 1 
VNT EIS EE I, (17 
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we find by the Thomas (1937) variational treatment the following two relations for th. 
constants y, g and и: 


u 


y — (ки — 3M) + (ха + Бхуи + и)", 
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wherein zy denotes the magnitude of the bound neutron wave vector x м: the latter is given by 
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| 2.518MeV 2MMr 
AN = z 


2x 2.518 MeV \ * 
h2 M+ Мт 


A | а | 
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The postulate of charge independence of nuclear forces suggests the assumption of 
equality of the neutron-target nucleus and proton-target nucleus interactions. Hence, im 
our calculations we take the potential Vp, to be of form (17). Using the momentum spaee 
representation for the functions Vpy, Vpr: Balp — ту). xT — Fr) we obtain, by the Dalitz 
(1951) and Lewis (1956) formulas, the following expressions for the matrix elements: 
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where f, and 7,A, are given as follows: 
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o = [С 
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he formulas for Е, and $54, are analogous to those for £4 and 7,A,; however, the symbol а; 
therein should be replaced by f. Tables I and II give the results of calculations and the 
experimental data. 


TABLE I 
| Calculated differential cross section | 
| 2 + 1 do | barn Experim. diff. 
| н | | | ST 0D pe cross section 
| ч A TE. Pur A NI те d barn 
ange o tentials VNT, En] rox 
| of the s entia NT: V pr ол 1.540: ater 
210-13 ст. | 3x1079em. | 47x10?em. | бх 10719 cm. 
0—09 0.80 | 1.02 1.46 | 1.55 0.04. 
a 10° 0.70 | 0.86 1.20 | 525 | 0.014 
0==9205 0.50 | 0.56 | 0.71 | 0.71 0.001 
9=40° 0.18 | 0.17 0.17 | 0.16 0.004. 


Table I contains the angular distributions of protons from the reaction 'N(d, p)'?N*. 
Here, the symbol !9N* denotes the excited state of the N nucleus, which is formed by 
the capture of a neutron in the 25 level. These angular distributions correspond to four 
ranges г = l/u of the potentials Vpr and Vyr. For each range of the potentials Vpr and Vyr 
the four values of the stripping differential cross section fix the shape of the angular protons 
distribution. These four values of the differential cross section are given for angles 0 = 0°, 10°, 
20°, 40° which the vector E, subtends with the vector b, In the last column of table I. the 
experimental values of the de differential cross section are included. From Table I 
it is readily seen that the experimental angular proton distribution is much more strongly 
peaked in the forward direction than the theoretical one. The same remains true even if 
we disregard the Vpr interaction. 


TABLE П 
Vli 
Values of the ratio ef prli > 
Range of the potentials (f |Vpnli> 
Урт, PNT Angle Angle Angle Angle 
6=0° 6=10° 0=20° 0=40° 
2х10 13 cm. 0.55 0.59 0.70 1.16 
8x10 !* em. 0.59 0.64. 0.76 0.29 
47x 10712 cm, 0.81 0.84 1.01 1.64 
6x 1071? cm. 0.81 0.86 1.04 1.68 
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From the Table П it can be concluded that, in the stripping reaction, the proton- target: 
nucleus forces are at least just as important as the neutron-target nucleus forces. The Born: 
first approximation tells us even more. Namely, for smaller binding energies of the captured: 
neutrons (e.g., for ey = 1.26 MeV) the matrix element < f |V pr|i> in eq. (10) can be quite 
decisive in fixing the shape of the angular distributions of the emerging protons. This situation 
arises on account of the resonances which occur at this small binding energy, and which 
sometimes yield monstrous stripping differential cross sections. The fact that no such effect 
occurs experimentally can serve as a hint that the capture process takes place on the edge 
of the nucleus only. 
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DIPOLE ARRANGEMENT IN PEROVSKITE-TYPE 
) FERROELECTRICS 


By A. JASKIEWICZ 
Institute of Experimental Physics, Wroclaw University, Wroclaw 
AnD H. KONWENT 
Institute of Theoretical Physics, Wrocław University, Wrocław 
(Received September 7, 1960) 


A method of investigating the dipole arrangement when perovskite-type crystals trans- 
form from non-polar to polar state is proposed. Thermal fluctuations cause thermal dipoles to 
be formed when the crystal begins to undergo transformation. A thermal dipole thus formed 
gives rise to an electric field in its surrounding favouring some directions for displacement 
of neighbouring ions. Discussion of ABO, substances with ferroelectrically active B ion indicate 
that there is one possible ferroelectric arrangement and three antiferroelectric ones. The ferro- 
or antiferroelectric arrangement is shown to depend on the polarizability of the А ion and on 


the lattice constant. 


1. Introduction 


'The structure of perovskite attracts much attention owing to its simplicity. The chem- 
ical formula of substances with this structure is 4BO 3, where A is a di- or monovalent 
metal and В a tetra- or pentavalent one. It is assumed that the ionic binding is more pro- 
nounced than the covalent one. Among these substances, there are ferroelectrics as well as 
antiferroelectrics. 

A ferroelectric crystal is characterized by the presence of spontaneous polarization, i.e. 
all elementary dipoles are aligned in one direction. An antiferroelectric crystal exhibits no 
spontaneous polarization. In an antiferroelectric, dipoles oriented in the same direction form 
chains along the direction of polarization, but two neighbouring chains are polarized anti- 
parallel. Such chain polarization, results in the total absence of macroscopic polarization. 

It is commonly assumed, that the thermal dipoles in ABO, crystals are formed by 
ion displacement. When all the ions are displaced in the same direction, thermal parallel 
dipoles are formed. If neighbouring ions are displaced in opposite directions, antiparallel 
dipoles are formed. An ion leading to the formation of a dipole is termed ferroelectrically 
active. In the materials under consideration, an A or B ion can be ferroelectrically active 
(281) 
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(Venevcev and Zhdanov 1956). Ion displacement takes place along one of the axes, parallel: 
to a cube edge, a face diagonal, or a body diagonal. 

The aim of the present paper is to investigate the mechanism of formation of ferro- 
and antiferroelectric arrangement in crystals of the perovskite-type structure, when only: 
the В ion is ferroelectrically active. 

Assume the crystal under investigation to trasform from the non-polar to the polar: 
phase. At sufficiently low temperature, a B ion is displaced and thus a thermal dipole is: 
formed. This dipole produces an electric field in its neighbourhood. The electric field in-- 
duces elastic dipoles in nieghbouring А and O ions. The electric field produced by the elastic : 
dipoles and by the thermal dipole acts on the nearest B ions and so facilitates the formation ı 
of further thermal dipoles. The direction of the B ion displacements depends on the directioni 
of the electric field and therefore the orientation of the thermal dipoles formed depends oni 
the direction of the electric field. By evaluating the electric field at the В ion sites, we cani 
obtain the thermal dipole arrangement and the state of the polar substance. 


2. Polarization of the elementary cell 


In Mason and Matthias’ theory (Mason and Matthias 1948) and in the theory of Piekara 
(Piekara 1954), the B ion is displaced into one of the six positions of minimum potential 
energy lying on the mutually perpendicular axes. In the cubic phase, the B ion occupies the 
symmetric position in an elementary cell, because its thermal energy is greater than the 
potential barrier U between the minimums of potential energy. 

At temperature T~U/k, owing to thermal fluctuations, one B ion may be displaced to 
ап asymmetric position, giving rise to an electric dipole. This dipole may induce elastic 
dipoles in the surrounding ions of the elementary cell. Such a polarized cell acts on the 
neighbouring ones as a polarizing agent. 

We consider an elementary cell of a perovskite-type single crystal with cubic symmetry, 
In the lattice, the A ions occupy the corners of the cube, the B ions are at the cube centre, 
and the O ions are at the face centre (i.e. we have a simple cubic lattice of 4 ions, face centred 
by O ions, and body centred by В ions, Fig. I). 

| 

It is assumed that the В ion is displaced along the z-axis in the positive direction, 

giving rise to a thermal electric dipole with moment m =m, Throughout the present paper, 


Fig. 1. Elementary cell of perovskite-type cubic crystal 
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this dipole will be treated as a point dipole. In evaluating its electric field, we make use 
of the well known formula for the field of a point dipole 


+ Sr т) — mr 
р S) 


5 (1) 


4 
The induced moment p; in the ion of type i at the point (x, у, z) is computed from the 

equation 
p (x, y, z) = a; Em, y, г), (2) 


where œ; is the polarizability of an ion of type 7. The coordinate system is shown in Fig. 1. 
The coordinates of the ions in an elementary cell with lattice constant a are given as 


follows: 


A3 
a a a a а w 
= : 9 
АА 5 2 Г 48 2" 2.5). (3) 
a а 
01 (o 0, 3 04 (- p 0, o) 
a a 
оз (оо, - $). os (0.4.0), 
а а 
os (4.0.0). os (o. - 5.9. 
B (0, 0, 0). 
The fields acting on the ions of type A and О are 
Вт. 
А 43 A4 A5 A6 "47 jo reca 
үй EX pua qt pe Ez E; E£ E? T 3y8 s 
AY 42 _ p43 44 _ pA5 р46 КАТ dan 8m; 
ali; e Е Е, ^y y y y 3y3 3" 
EA! sr Е42 UE Баз = pa pe = pee oe ps 0. 
ЕФ is p? e E93 P EB pU p 0, (4) 


a Оз O4 O5 O6 . 
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Or. роз _ 10т. 
E; SU IE 
8m, 
ОЗ _ p]pO4 _ 1057 = 106 — kesh, = 
Е, E Е, къ [p ДЕ аЗ 


Denoting the ionic polarizabilities & ,. Xo for the ions A and О, respectively, we obtain 


the following values for the dipole moments induced in the ions of the cell under considera- 


tion: 
A5 A6 A7 48 Sma 
AL EE AS 7 а. Aa Ар oe аа. an = 
Рх = Px > Px Рх Рх Рх Px Px 43/303 
m, 
A A2 A3 44 АБУ AG m qm AS z 
ру pui = Pon = Poem р нр: p ПА 


~ 74375 a8” 


O1 — р02 _ 03 — „O4 O _ 06 0, (5) 


NE 0 l6m. 
Pe = Pg = по = 
x m. 
: : j: 
Pe pei роб ро — д0 PE 


where рї is the х-сотропепі of the dipole moment induced in ion 41. 


3. Electric field at B ion sites 


In order to evaluate the polarization induced by the polarized cell in the surrounding 
ones, we take into account the interaction between the nearest neighbours, The polarized cell 
is surrounded by six adjoining cells with the faces in common. Polarization of each joint 
face is already known, the latter being part of the polarized cell. The cells are numbered 
as follows: the cell adjoining to the face centered by the Ol ion gets number one, that adjoining - 


to the face centered by the O2 ion gets number two, etc... The remaining cells have polarized 


corners and edges. The influence of these corners and edges on the cells they have in common 


is relatively small and therefore in further considerations may be neglected. 
The polarization induced in the surrounding cells by the polarized cell may be easily com- 
puted if the electric fields are known. Therefore, now proceed to compute the electric fields 


acting on the B ions at the symmetry positions in the cells numbered 1, 2, 3, 4, 5, and 6 due 


to the dipoles in the polarized cell. With regard to the symmetry conditions, we compute 


separately the electric field acting on the cells 1 and 2 and on the remaining ones. 
The electric field strength at the centre of cells 1 and 2 due to the dipoles induced in 


ion A of the zero cell (polarized cell) is given by 


1 | е 912 хат, 
Ew P о IE: | (6) 
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The electric field strength due to the dipoles induced in ion О at the centre of cells 1 and 
2 is given by 


EB, = Ey =0, Б,- 2562975, (7) 


The electric field at the centre of cells 1 and 2 due to the thermal dipole of the zero cell is 


Е | и Mm 
E,=Ey=0, B=. (8) 


а 


(8) 81 


(e) 82 


Fig. 2. Configuration of the B ions nearest to the primarily displaced В ion 


The total electric field at the centre of cells 1 and 2 is 


O о? а та 


а Ж О PM (9) 


We are now able to compute the electric field at the centre of cells 5, 4, 5, and 6 given 
rise to by the dipoles induced in the O ions: 


Due Жесс? (10) 


a 
The field due to the dipoles induced in the A ions is 
E, = Е, = E, =0. ~ (11) 


'The electric field due to the thermal dipole of the zero cell is 
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The total electric field at the center of cells 3, 4, 5 and 6 is 


( m 
EP ДОБО, ee = (6:20 A As (13) | 

We see that, at the sites of the B ions in cells 1 and 2 electric fields are acting which ı 
tend to displace the B ion in the positive direetion along the z-axis, and others tending to | 
displace it in the negative direction. The field acting in the positive direction is produced by : 
the dipoles induced in the O ions (eq. (6)). Hence the resulting field depends on the same : 
parameters as eqs. (6), (7) and (8). The Bl ion is displaced in the positive direction if the 
following condition is fulfilled: 


E® > 0. (14) 
For the displacement of the Bl ion in the negative direction we have 
ЕФ <0. (15) 


Considering the electric fields acting on the B ions in the side cells, we see that the 
field due to the dipole induced in the O ion as given by eq. (10) is oriented in the positive 
direction along the z-axis, whereas the field produced by the thermal dipole of the zero cell, 
as given by eq. (12), is oriented in the negative direction. The electric field produced by 
the A ions is zero, as is obvious from eq. (11). The electric fields given by eqs. (10) and (12), 
respectively, compete to displace the B ion in the positive or negative direction. The B ion 
will undergo displacement in the positive direction if the following condition is fulfilled: 


ЕФ > 0. (16) 
The B ion is displaced in the negative direction if 
ЕФ = 0. (17) 


The four conditions (14), (15), (16) and (17) yield four possible arrangements of the 
dipoles in the ferroelectric and antiferroelectric state. For a given crystal structure, we obtain 
the ferroelectrie state (Fig. 3a) if the conditions (14) and (16) are fulfilled. Pairwise, the‘ 


Fig. 3a Fig. 3b Fig. 3c Fig. 3d 
Fig. 3. Possible thermal dipole configurations for perovskite erystal when B ion is displaced from its original 
position 
Fig. За. EY 50, ЕФ у Fig. ЗЬ. EY 50, ЕФ) <0, 


x x 
Fig. 3c. Е) = EC? >0, Fig. 3d. EM <0, Ee <0. 
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conditions: (14), (17); (15), (16); (15), (17) yield the antiferroelectric arrangement as pre- 
sented in Figs. 3b, 3c and 3d, respectively. The three configurations 3a, 35 and 3d were 
discussed by Kinase (Kinase 1955) using the local field method. 

The resulting fields EP and ЕО will cause permanent ion displacement if the energy 
of the thermal dipole thus formed in the electric field EC? or E is of the order of the thermal 
energy. 


4. Applications 


For ABO, substances, the polarizability of the O ion is хо = 2.38 x 10-4 cm? (cf. Slater 
1951). Since the polarizability of the oxygen ion is known, we can write the conditions for 
a ferroelectric arrangement of the crystal in the form dependent on the lattice constant and 
the polarizability of the А ion. The conditions, in this new form, are readily applied to any 
particular ABO, substance. 
| From eqs. (14) and (16), using eqs. (9) and (13), we easily derive the inequalities 


which must be fulfilled for the substance to be in the ferroelectric state: 


ад « 48.4 x LO cm? (18) 
a SSK lO cm (19) 


In BaTiO, and KNbO, crystals, the B ion is known to be ferroelectrically active; there- 
fore, our results may be applied to these crystals. Since, in these crystals, the polarizability 
of the А ion and lattice constant fulfill conditions (18) and (19), these substances are ferro- 
electric. 

An antiferroelectric arrangement results if, for a given crystal, inequality (18) or (19) 
is not fulfilled. As the polarizabilities of all known ions are of the order of that of oxygen 
ion, it is highly likely that condition (18) should be fulfilled. Therefore it may be anticipated 
that if the thermal dipoles are formed by the В ions only, ABO; crystals will present the 


ferroelectric state. 


5. Conclusion 


In this paper, we proposed a method of investigating the dipole arrangement occurring 
when a crystal transforms from the non-polar to the polar state. We discussed the dipole 
arrangement which will develop then. As transformations to a polar state set in, one ion is 
displaced from its original position forming a dipole and giving rise to an electric field in 
its neighbourhood. 

From our considerations it may be inferred that the characteristic features of ordering 
of thermal dipoles depend essentially on the elastic dipoles induced in the A and О ions. 
Therefore, the ferro- or antiferroelectric arrangement will depend on the configurations 
and polarizabilities of the ions. The induced dipoles will not, however, affect the total electric 
moment of the elementary cell, since its resulting elastic electric moment is zero, as is easily 


inferred from formulas (5). 


288 


It is shown that, in crystals of perovskite-type structure, there are four possible kinds 
of the thermal dipole arrangements, for thermal dipoles formed by the B ions. One arrange: 
ment is a ferroelectric one and the three others are antiferroelectric. 

In polar ABO, substances, with В ion ferroelectrically active, the ferroelectric arrange- 
ment can occur if the polarizability of the 4 ion is smaller than 48.4 x 10-24 cm3 and if simul- 
taneously, the lattice constant is smaller than 5.35 x 10-8 em. In known substances, the condi- 
tion for the ferroelectric arrangement is fulfilled. It is expected that further ferroelectric 
and antiferroelectric arrangement will be found when considering ABO, substances with 


the A ion ferroelectrically active. 
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TRANSPORT OF CURRENT CARRIERS IN n-TYPE INDIUM 
ANTIMONIDE AT LOW TEMPERATURES 


By J. KOŁODZIEJCZAK 
Institute of Physies, Polish Academy of Sciences, Warsaw 


(Received September 29, 1960; translated paper received November 25, 1960) 


The electric conductivity, Hall effect, and electron mobility are computed within the 
temperature range of 0° K to 100° K, account being taken of the non-parabolic structure of the 
conduction band, the ionic mechanism of scattering, and a weak magnetic field. The depend- 
ence of the effective mass on the electron concentration is derived. The transport equation 15 


solved by Me Clure's general method. 


1. Structure of the Conduction Band 


According to results due to Kane (1957), strong interaction between the conduction 
band and the valence band of light holes perturbs the parabolic character of either. Interaction 
with all remaining bands is considerably weaker and contributes anisotropic corrections to 
the energy that are negligible at 4? «€ 5 x 10-4 atomic units. The absence of anisotropy in the 
cyclotron effective mass seems to corroborate the spherical band model (Dresselhaus et 
al. 1955). Measurements of absorption on free carriers (Spitzer et al. 1957) point to an effective 


of the electrons depending on the concentration of the latter. An identical result was 


mass 
e experimental 


obtained from measurements of the Faraday effect (Smith et al. 1959). Th 
dependence of the effective mass on the electron concentration confirmed the results of 


Kane's theory. The various energy bands are described by the equation 


ele e (eet Д) Kk? P(e’ + 84) = 0 (1:1) 
wherein 
\ пада 
d а 2mo 


at k=0, and A denotes spin-orbit splitting. The constant accounting 


ё, is the energy gap | 
nd is defined as follows: 


for coupling between the conduction band and valence ba 
7 
Duc A ГАР: (1.2) 
та 


wherein S and Z are the wave-functions corresponding to either band. For A>kP,e,, the 


(289) 
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expression describing the conduction band assumes the form 


8132. ] 8 ” 
t= i : = ie (4+ В чә) | ; (1.3 


Eq. (1.3) yields an accurate description of the structure of the conduction band at 2 <5 > 
x 10-4 atomic units, corresponding to electron concentrations of less than 1019 cm 3, 


2. Transport Equation 


Assuming the existence of a relaxation time and the presence of an external magnetir 
field, the equation of transport can be written in the form 


tle i) grad, f+ I —0 (21 


with the following notation: 


> > 


E — electric field strength, H — magnetic field strength, v — electron velocity, f — distri: 
bution function, f, — distribution function at thermodynamic equilibrium, т — time 0: 


relaxation. 
The transport equation (2.1) was solved according to Mc Clure's method, yielding tha 
following components of the conductivity tensor: 


2e? 2 ы | 


wherein S are Ше components of a tensor, given by 


۹ у, 2|v.(m)|? 
1+ (тот)? 


m=1 


3 + umet 
1+ (тот)? 1+ (тот)? 


SE Е x(m) vy(— т) + г, (- m) (т). vx(— m) vy(m) — vx(m) vy(— 2 
m=1 

nu | 
Sse = 280) + М I чы (2.3) 


m=1 


w is the cyclotron frequency, and v,(m) — the coefficients of the Fourier expansion of the: 
velocity. Sya Szy and Sy, are obtained by varying the indices of v(m). The Fourier coefficients 


v(m) are b by er the equation 
Germ | 
hee = т Dos dal (2.4) 


The velocity of the electron is given by 


ПЕРСЕ баща 8 E 
ae тена (4+ er) р}, { (2.5) 
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With the notation 
he A 8 с сад 
h2 E БЕ 3 (4: 3 вр) r] =m* А (2.6) 


. (2.5) сап be rewritten аз follows: 


U= x 


(2.7) 


m 


ence, the electron can be considered to be free if its mass is assumed to be given by eq. (2.6). 
The solutions of eq. (2.4) are the following: 


e ; : 
жыл Ба, райе“. iot —iot 
Ux aH 9 9/1 + е 
EN wat ; 
ос но. kp [et — e~] са 
vith о denoting the cyclotron frequency as given by 
eH fh? 46, Вер} pil = еН 
ups 43 3 (в Sup B am (2.9) 
p is the projection of the wave vector on the plane perpendicular to the direction of the 
nagnetic. field. The following relations hold: 
Е, = (Б) k= (k + ke)? (2.10) 
3y (2.8), the following Fourier coefficients are obtained: 
СПК! 
ы) =ы(—1) = 5 0 > hy 
"LL ü 1 
vy(1) = —v(—1) = aH ® 5 kp (2.11) 
h 
v; (0) = = Он 


d] other coefficients vanish. 
With eqs. (2.2) and (2.3), the components of the conductivity tensor can now be written 


s follows: 
E Uy = C i} ر‎ ris pee dk 


бу ME E c tdk (2.12) 


0, =0= ДЕ? 
0 


ае о асан 
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wherein 


ма fae а RE 
127 де m* 


Subsequently, only the case of so-called weak magnetic fields will be considered, a 


determined by the condition то 1. The conductivity tensor components are now given bl 


= 
O xx = Oyy Nr c | аак (2.82 
0 
x А 
Oye = — o, = C| Атомна! 
0 


е 2 


3. The Conductivity Tensor at Low Temperatures 


The scattering mechanism predominant throughout the temperature range of 0—100° F 
is that of scattering on ionized impurities. The relaxation time for the case under consider 
ation is (Barrie 1956): 


hz? de | x i Ge 
= ees, г ee 
VS Grete die | Becas 
Е = 2 (3.1 
F 


with N denoting the concentration of the scattering centres 


х —the dielectric constant, and 
q —the screening constant. 


The scattering potential is given by 8 


V(r) вера ее (3.2 
xr 

Eq. (3.1) is of a higher degree of accuracy than the Conwell-Weisskopf formula, as it takes 
into account the screening of the scattering centres by the electrons. The sc 
is a function of the electron concentration. In order to facilitate subsequent computations. 
it is of advantage to proceed to the integration of the expressions of eq. (2.13) 
‚ Assuming zero energy at the bottom of the conductivity band, eq 
as follows: 


reening constant а 


over the energy. 
‚ (3.1) can be rewritte 


mm 
fo 


5 h2%2 1 £a oan 
ё = 2m, + 5 €g + 3^ f. — е (3.3) 


The constant P can be determined by comparing the effective cyclotron mass (2.6) wit 


the value measured at cyclotron resonance. Such measurements carried out with samples 


of a very high degree of purity and at low temperatures (Dresselhaus et al. 1955) yielded 
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еспуе mass of m, = 0.013 mo. Considering the high degree of purity of the samples and 
1e low temperatures, this can be assumed to be the mass of the electrons at the bottom of 
1e band of conductivity. 


Putting Ё = О in the effective mass of eq. (2.6), and equating m* and m,. we have 


S Mg Un 


22 Seen 
quieta mont, d 
olving eq. (3.3) in Ё vields 
m = yy, 
k n mt 2yu —V1+ Ax(u — u?) }” (3.5) 
'herein 
€ Mn 
===—: L= 
: Eg А Mo 


1,15 the cyclotron mass of the electrons at the bottom of the conductivity band, i.e. at & = 0. 

Considering that и is a small quantity amounting to 0.013 (the latter results from meas- 
wements of m, at cyclotron resonance), the root in the brackets in eq. (3.5) can be expanded 
n a power series. This procedure is fully justified, as the inequality 


LUI. 

& == PASS 3.6 

4u — u?) с, 

s fulfilled even in samples containing the greatest amount of impurities. Indeed, at e, = 

= 0.23 eV the right hand side of the inequality is of the order of 6 eV, whilst the left hand 
ide, at an electron concentration of 108 em-?, amounts to 0.16 eV. 


Expansion in series yields, in the first approximation, 


fies ee ХО + x)]* (3.7) 
vhence 
A dk = a [x(1 4-30] ^ (1+ 2x) dx (3.8) 
Phe effective mass of eq. (2.6) now assumes the form , 
m* = m (L+ 2%) (3.9) 


Phe relaxation time of eq. (3.1) assumes that of 


x? (2m,)^^ [e(l + y)] ^ 


3.10 
aNetB 1+2% ыл 


т = 


vith 


8m, , 
B= 0+8: E= aa ey (3.11) 
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In computing the screening constant g, the method proposed by Dingle (1955) is utilized: 
however, it was necessary to make allowance therein for the non-parabolic structure о 
the conductivity band. 

The concentration of the electrons is given as follows: 


Qn [ast (3.12: 
by (3.17), this yields 
4 27 Ws jt °/ 2 3/, 
E | "ur | LOB, т (3.18 


The (В, n) are integrals defined as follows: 


= | | 3 (x + fr" (1+ 283) dx 
dx, 


mx RT 


with ky denoting the Boltzmann constant, and © — the Fermi level. 
The electron concentration about a centre of positive charge is 


4 алт TV 2 ар, ms | 
P a | M Jure d (3.14 
with 
; y e 
ye = pla qc a 
kt ur“ 


The excess space charge density is given by the following expression: 


4 2 „КТ 3 t 
م‎ еа = — 0 5 | ns Va (LOB, + V*) — 058, эр} (3.15) 


The potential energy of the electron in the Coulomb field of the centre can be assumed tc 
be small as compared to the kinetic energy. This is justified inasmuch as the probability 
for an electron passing at a very small distance from the centre is small. With this assumption. 
the quantity L(*(8, n+ V*) can be expanded in a series, thus 


4 x T a 
LEB. 3 у) I. л) + [S Kr (3.16) 


aV* | 


The density of space charge is now 


Are? (2m, !: 
0= : = - (FT) ] "a " (D. n) d (3.17) 
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Ра 


oisson’s equation 
4л 
Ру = — — 
та: (3.13) 
ust be fulfilled. Substituting the scattering potential of eq. (3.2) and оп making r tend 
o infinity, we have 


a Поле? л гуя 


: "E = LY B, m) (3.19) 
Eq. (3.11) now assumes the form 
xhe(l-+ y) 
em, T)^ LE 120) 


Because of the weak variability of the coefficient B appearing in the formula of the relaxation 
time, it can be placed before the integral with an energy value corresponding to the maximum 
sub-integral expression, 
5 
х«Ь(ЪТ)% Ld* 
e*(2m,)* Li? Le’ 


Е m 
Smax 


(3.21) 


On ground of the foregoing considerations, the components of the conductivity tensor (2.13) 
assume the form 


Oxy = 0,, —0,, — 0 = GL, (3.22) 
gy, = — Oy = GV y Ll 
wherein 
16 е 
G= = л pu та lly 

a 2x (Kk T)? 

9 zNe*B 
e2H? 


Ayr 
deni 2т. КТ 
The quantity В is defined by eq. (3.11). 


4. Electric Conductivity and. Hall Effect 


From the considerations of the preceding section, the electric conductivity апа Hall 
effect in a transversal magnetic field can be computed. 'The current density depends on the 
electric field strength as follows: 


J, — o, E, 4 o, E, (4.1) 


xx H 
Jy = ДР uy gy Es 
Putting dy = 0, and with respect to с, = — 0,45 0y, = Oxy, We have 


ERIE J, = oq; (4.2) 
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R is Hall’s constant, and бн — the electric conductivity in a transversal magnetic field. 


Both quantities are given by the conductivity tensor components: 


2 2 
) Oxx T Oyy - 
p н ы (4.3) 


2 a? 
H Oxx T Ovx Oxx 


R 


For a vanishing magnetie field we have 


зул Ба па 


сър = 4.4) 
BR, 8e h? иу ( ) 
or 
ПРА D 1 
Ro em TENE (4.5) | 
and 
я = GU. (4.6). 
Aen e Ва 


The numerical values of the coefficient = ў corresponding to different В are given 


3 9 

(LI 3)? 

in Fig. 1. It is noteworthy that, if —0 is put in the foregoing formulas. the respective 
expressions for the parabolic band are obtained: 
3 1 ЕЁ, 
4 еп (KF)? * 


The functions F, in the foregoing formulas are Fermi integrals defined 


Ко = с = 3GF, (4.1) 


as follows: 


BRÜHL, 


Fig. 1. Dependence of the coefficient rY Q3. on n, for different values of DOM 
ша 


ге 
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= 
= . » 
e numerical values of L; have been assembled in Tables I, П and IM. 
TABLE I 
со 
d ж + 8; „2\3 
Values of D E | ү тре) 
се dx | (1 + 2px)? 
0 
8 с - = 
0 0.05 0.10 0715 0.20 | O25 0.30 
= rer 2. = | 
0.109 0.097 0.094 0.094. 0.0966 0.0995 0.103 
0.798 0.710 0.687 0.689 0.704 0.725 0.751 
2:115 1.881 1.822 | 1.829 1.868 1.927 1.997 
5.409 4.793 4.645 4.669 4.777 4.935 5.124. 
12.94 11.50 11.16 11.25 ПШ 55 11.98 12.47 
28.33 | 23:31 24.63 24.96 25:476 26.83 28.06 
103.6 89.56 88.54. 91.48 | 96.19 101.8 108.0 
TABLE II 
се 
9 К Af | ст Вх? NE 
Values of p = Дона: dx 
С dx | 0+ 202) 
0 
? p «t З я 
0 0.05 0.10 0.15 0.20 0.25 0.30 
4 = 
| 
E 0.957 0.516 0.378 0.313 0277 0.255 0.240 
وو‎ 7.042 3.790 2.779 2.304 2.038 1.873 1.765 
=] 18.96 10.22 1.491 6.211 5.496 | 5:0555 4.761 
0 50.33 26.81 19.66 16.29 14.41 13.26 12.51 
+] 129.7 68.51 49.97 41.39 36.66 33.17 31.91 
+2 318.4 163.8 119.0 98.66 87.57 80.89 76.62 
+4 1560 746.7 535.9 447.0 400.7 374.1 350-4 
TABLE Ш 
со 
3 2 df E 
Values of iy 3 | | (x + Вх2) ldx 
0 
к 0 0.05 0.10 0.15 0.20 0.25 0.30 
7 ; 
—4 0.0161 0.0194. 0.0228 0.0265 0.0303 0.0343 0.0385 
== OMS 0.138 0.163 0.184. 0.217 0.246 0.277 
=I 0.291 0.345 0.408 0.476 0.548 0.623 0.702 
0 0.676 0.833 0.996 1.170 1.354 1.548 Gy 
+1 1.396 1.750 DEO 2.534 2.965 3.420 3.889 
+2 2.502 3.233 4.032 4.893 5.812 6.787 7.813 
- БО ДА 8.198 10.84 13.73 16.85 20.17 23.70 


298 


5. Strong Degeneracy 


The non-parabolie structure of the conductivity band can produce measurable effects 
in the case of samples with a high amount of impurities only. High concentrations of the 
latter, especially at low temperatures, cause strong degeneracy of the electron gas. Indium 
antimonide, which presents a small effective mass of the electrons, is а material wherein 
degeneracy can be especially strong. The Hall constant as measured by Putley (1959) 
in samples of different concentrations of the impurities is temperature-independent through- 
out the range of 2° K to 100° K. The absence of activation energy resulting can be explained 
by the existence of an admixture band overlapping the conduction band. The presence 
of electrons in the conductivity band warrants strong degeneracy even at the lowest tempera- 
tures. Obviously, the degree of degeneracy is dependent on the concentration and tempera- 
ture, and can be determined by solving the following equation with respect to 7: 


2 LE = 1.24 x 10-135 7l (5.1) 


n is the reduced Fermi level. 
"Strong" degeneracy is meant to denote the case of 7 > 4. At 2°K, ай samples with con-. 
centrations exceeding 10! cm? already exhibit strong degeneracy. It can be proved that, . 


d : 
Karimas (ee P. > ô(x — 7). Introduction of the delta function makes the exact соте: 
К 


putation of the integrals ZA possible. Eq. (3.13) yields 


¢=£(V4—-1) (5.2) 


wherein 


па (5.3) | 
NTL em, 
On passing to infinity with с, and expanding the root in (5.2) in a power series, the respective | 
expression for the parabolie band is obtained: 


== Жа МА 
ct] 8m, ‘ (5.4) 


The electric conductivity of eq. (4.6) now assumes the form 


та Р 
1 3 һЗх?п? 
а = 


2 \ x) ети (5.5) 


ге] effective аве RFA E үнү = ; 
Herein, the effective mass m* is given by the concentration: 


E МА 2 м 
пе ез та | 1 - x ed ns (5.6) 


л Eg Mn 


The quantity B is determined by eq. (3.11), except that now 


8 


л Дет? 


Sr 


(5.7) 
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If Analogical procedure with respect to the Hall constant of eq. (4.5) yields 


Ry = — — 5.8 

TE (5.8) 

The dependence of the effective mass (5.6) on the concentration is shown in Fig. 2. The 
computations were carried out for e, = 0.23 eV and m, = 0.013 mọ. It will be seen from 
the graph that, as the concentration exceeds 1016 cm-3, the effective mass increases steeply. 
Experimental values of the mass from measurements of the absorption on free carriers 


E 0.00277 0.0124 Frergy 00496 0.164 ev 
5 
4 
m* 
m, 
3 
2 
; ا ا‎ 
107 


10% 10" 


104 1015 1016 м 


Concentration cm 


Fig. 2. Effective mass m* versus electron concentration. 


(Spitzer and Fan 1957) and of the Faraday effect (Smith et al. 1959) are also shown. Compari- 
son of the conductivity of eq. (5.5) with the experimental results is less conclusive, since 
this involves an unknown parameter, namely, the concentration of scattering centres JV. 
Assuming the electrons to be scattered on ionized donors and acceptors, we have 


N = ng p, (5.9) 
The condition of electric neutrality of the crystal is fulfilled: 
пр =p, (5.10) 


with n denoting the electron concentration in the conductivity band, 

p — the concentration of holes in the valence band, 

ng — the concentration of ionized donors, and 

р, — that of ionized acceptors. M 
The assumption of p —0 at low temperatures is reasonable, as all acceptors are ionized by 
electrons from the donors. Hence, the concentration N amounts. to 


N = 2p, (5.11) 
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Finally. we obtain the following expressions for the electric conductivity and electron mo- 
bility: 


1/3) Ber? 
TE eB 5.12 
М 18 | л | (N; + n)e?Bm*? ( ) 

3 \* ЛЗх?п 
ЕУ | 3 (N; + n)e3Bm*? (5.13) 


wherein: 
e. 
N; 8 Pa 


Graphs of (5.12) and (5.13) are shown in Figs. 3 and 4. The computations were carried 
out for different values of N, at 4 —17.5. The experimental points are those derived from 
the references. The electron concentration for each sample was determined from eq. (5.8) 
(from the Hall constant as measured at low temperatures). 


$27 ст"! 


Conductivity 6 


102 


Fae AC 
ar. m i 
ДИ | 


105 70% 


1017 10 8 1019 


1020 
Concentration n 


cm 
Fig. 3. Electric conductivity versus electron concentration, for different values of N;. 
lhe experimental points are taken from the papers by Putley (1959), Taue and Matyas (1955) and Frederikse 
and Hosler (1957). 


In samples of a high degree of impurity, the electric conductivity and mobility are 
temperature-independent throughout the entire interval of impurity conduction. In samples 
that are of a somewhat higher degree of purity, the conductivity and mobility are found to 
increase as the temperature is raised. This latter effect is due to the vanishing of degeneracy. 
Ihe temperature dependence of the conductivity is described by eq. (4.6). Fig. 5 brings 
graphs computed from that formula for three of Putley’s samples. The graphs corresponding 
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mobility д 


10* 


103 
Concentration n 


Fig. 4. Mobility versus electron concentration, for different values of N.. 
The experimental points are those for conductivity. 
The dashed line corresponds to the parabolic band. 


to samples of a rather high degree of impurity are in good agreement with the experimental 
data. The results for the purest sample C48/4 are less consistent, although the shape of the 
graph is conserved. The graphs were computed on the assumption that the electron concen- 
tration is temperature-independent for a given sample throughout the range of impurity 
conduction. In determining the concentration of the electrons from Hall's constant as 
measured at the lowest temperatures, one can compute the temperature dependence of the 
Fermi level from eq. (5.1) and then proceed to compute the graphs of Fig. 5 from eq. (4.6). 
Lack of consistency with the experimental data in the case of the samples of highest purity 
can result from an incorrect evaluation of the concentration. The source of error can reside 
in the fact that, for the magnetic fields of 2000 Oe utilized by Putley, the classical theory 


of the Hall effect ceases to be valid, as then 


ho > 
with w denoting the cyclotron frequency and С — the Fermi level. The foregoing inequality, 
if fulfilled, sets a limit to the applicability of the classical theory. For a concentration of 
1014 em-3, Fermi's level is given by a value of 5.98 х10 7 eV, whereas the Landau splitting 
amounts to hw — 1.74 x 10-3 eV. The quantum effect becomes appreciable here. In the case 
of samples of higher carrier concentration, the quantum limitation sets in at higher values 
of the magnetic field. strength. 

The author wishes to thank Professor Dr L. Sosnowski for directing the present inves- 
tigation and for his valuable hints. The author is indebted to Dr Dr M. Suflezynski, J. Ginter 


Temperature °K 


Fig. 5. Electric conductivity versus temperature. The experimental points are from Putley’s paper. 


and J. Rauluszkiewiez for their critical remarks. The author desires to express his special 
acknowledgement for the effective help provided by Mr. B. Zaloga in the field of numerical 


computation, 
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LUMINESCENCE IN ELECTRIC FIELDS OF “DRY” OXIDE FILMS ON 
ALUMINIUM 
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М. JACHIMOWSKI* Ако R. DRAGON 


Institute of Physics, Teachers’ College, Opole 
(Received October 12, 1960) 


The electrical and luminescent properties of electrolytically formed “dry” aluminium 
oxide films were investigated. Cells consisting of an aluminium base with oxide layer and trans- 
parent semiconducting coating were subjected to a D.C. field. 

The light emission from the oxide layers was observed. 

The luminescence brightness and cell current were measured versus the applied voltage. 

The properties of dry aluminium oxide films strongly resemble, on the one hand, those 
of electrolytic cells and, on the other hand, those of typical electroluminescent cells. 

It is therefore assumed that in either case, viz., in that of dry cells and in that of electro- 
lytical cells, the light emission may be considered to represent a process of electroluminescence. 


1. Introduction 


The luminescence accompanying processes of anodic oxidation of Al, Ta, Zr and some 
other metals is a phenomenon which has been known and repeatedly described for over 
50 years (Güntherschulze 1906; Lavaux 1919; Dufford 1929; Forrest 1930; Sullivan et 
al. 1931; Güntherschulze et al. 1932, 1934; Betz 1935; Verwey 1935; Rummel 1936; Gu- 
miüski 1936; Anderson 1943). 

During the last decade and, particularly, їп the last few years, 


received special attention. 
Systematic, extensive investigations were carried out and a number of theories aimed 


at explaining the origins and mechanism of the phenomenon were put forward (van Geel 


et al. 1950, 1956, 1957; Dekker 1950; Dekker et al. 1950; Charlesbey 1953; Ruziewicz 1956; 
Smith 1959). 


this phenomenon 
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\luminium was the metal investigated most extensively. Aluminium plates or foils 

of various grades of purity, formed the anodes of a suitable electrolytic cell (usually contain 

ing an aqueous solution of boric acid or oxalie acid), the other electrode (the cathode) usu: 
ally being made of platinum. 

- "Let us recall some basic facts described in the literature, i.a. in the paper of van Gell. 
Pistorius and Bouma (1957), where they are presented in a concise and particularly clear 
way. 

| (а). The intensity L of the “anodic”? luminescence occuring during the formation of the 
oxide layer or in a completely formed layer depends on the actual thickness d of the layer 


and on the density / of the current flowing through the system: A1/ALO, electrolyt. 
L = al(exp Bd — 1), (1) 


х and Û being constants for a given sample. 
? 


(b). At a constant thickness d of the oxide layer, the current density / is a function: 


of the voltage V applied to the electrolytic cell. 


I = a,(espbV — 1), (2) 
a, and b being constants. 
By (1) and (2): 
Г 41, (3) 
L = a(exp bV — 1). (4) 


a,b and c being constants for a given specimen. 

(c). The system AI/Al,O,/ electrolyt is known to exhibit rectifying properties, the: 
blocking direction being that from (Al +) to (electrolyt —). At the opposite polarity, (Al—), 
larger currents flow through the cell. a gradual deterioration of the structure of the oxide: 
takes place, and the system loses both its rectifying and luminescent properties. A sudden: 
reversal of the polarity, from (Al +) to (AI—), produces short-lasting luminescence, termed | 
“the cathode flash”, on the aluminium electrode. The cathode flash occurs only when the: 
Al previously constituted the anode, and сап be observed even after ац interval of several! 
hours between the removal of the anodic polarization and its sudden reversal. 

A subsequent, sudden transition from the state (А1—) to the state (Al +). results ini 
а short-lasting “anode flash” superimposed on the normal, steady anodic luminescence. . 

The anode flash occurs when the Al previously constituted the cathode, provided the: 
currentless time interval separating the states (Al—), (Al+) does not exceed several minutes. . 

\n anode flash of low intensity also occurs when the anodic system is first short-circu- - 
ited and then submitted to the previous anodie polarization. | 

These facts, and many another described in the literature, suggest that we are dealing : 
here with a phenomenon of electroluminescence of the aluminium oxide (Ruziewiez 1956; : 
Curie 1953; Curie et al. 1954; Piper et al. 1958; Smith 1959). 


If this is so, aluminium oxide in the "dry" system: AI/Al,O,/ semiconductor or metal | 
can be expected to exhibit similar properties, / 
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This paper gives ап account of some experiments carried out on dry aluminium oxide 
lls. 

Fach cell consisted of an aluminium base plate covered with an electrolytically formed 
xide layer and a transparent film of “n” type semiconductor deposited on the oxide. 


2. Apparatus 


2.1. Luminescent cells 


Aluminium plates (of about 20 cm? area), generally of technical purity, were degreased, 
‚ashed in distilled water and oxidized anodically in 1—2% aqueous solutions of oxalic 
cid. Two methods of oxidation were adopted: the constant voltage method and the constant 
urrent method. 

In the first of these, a constant formation voltage is applied to the electrolytic cell 
ntil the continually decreasing current drops to the low and nearly constant value of the 
leakage" current. | 

In the second method, D.C. of constant density (in our experiments 2—5 mA/cm?) 
s passed through the system until the increasing voltage reaches the preselected value of 
ormation. 

Formation voltages ranging from 24 to 140 volts were used. The oxidized specimens 
ere washed thoroughly in running and distilled water and dried in air at room temperature. 


inally, part of the oxide area was coated with a transparent and conducting CdO or SnO, 


ilm. 

Samples intended for CdO coating were mainly oxidized by the constant current 
orocedure. 

Cadmium oxide films were deposited by the cathodic reactive sputtering technique 
n an apparatus similar to that described by J.S. Preston (1950). 

` Transparent layers of SnO, were fumed by the chemical deposition procedure (Fischer 
1954; Holland). Notwithstanding the temperatures required in this method (~500°C), по 
breakdowns were observed in the majority of specimens throughout the range of voltages 
ased in our experiments. The luminescent cell finally constituted a parallel plate condenser. 
Voltages were applied to the aluminium base and to the semiconducting film by means of 


an elastic spring contacting the latter. 
There appeared to be little difference in the behaviour of the cells prepared by either 


method. 


2.9. Luminescence 


Quantitative investigations of the luminescence of oxide films were carried out with 


ı F. Е. U. 19-M type photomultiplier in the simple arrangement of Fig. 1. 


The arrangement allowed for simultaneous measurements of both the intensity 1, 


of the photocurrent, which is proportional to the luminescence brightness L, and the inten- 
sity I, of the current flowing through the cell, as functions of the value and polarity of the 


voltage applied to the cell. 
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Aluminium plate 
| Oxide layer 
/ / 


Fig. 1. Block diagram of the apparatus. Sı, D. C. voltage source; 5,, supply for photomultiplier; V, voltmeter; 
C, voltage commutator; 4;, Az, microammeters; Е, photomultiplier. 


3. Experiments and results 


3.1. Anodic luminescence and flashes 


All experiments were carried out in air at room temperature. Preliminary visual obser- 
vations showed that the cells prepared as described in Section 2.1 emitted light under the 
influence of sufficiently high voltages applied in the blocking direction, (Al +). 

At the opposite polarity no steady luminescence was observed. Sudden transition from: 
the state (Al+) to the state (Al—) produces a cathode flash which, however, decays at 
a much slower rate than in the case of an electrolytic cell. A subsequent sudden return to the 
state (Al--) causes an anode flash, superimposed upon the steady anodic luminescence. 
This anode flash is also of longer duration than in the electrolyte. The anode and the cathode: 
flashes are mutually dependent in the same way as they are in the electrolyte. | 

Аз in the case of the electrolytic cell, the 


anode flash of low intensity can be observed. 
without previous cathodic polarization, i. e. 


when the system, upon anodic polarization, 
is first short-circuited and then subjected to the same polarity again. 


3.2. Relationships between Voltage, Current and Luminescence 


On applying various voltages V to the cells 


‚ cell-currents J, proportional to the current 
densities /, and photocurrents 7 


; were measured. The general trends of the relationships 
L), 14У) and L(V) were approximately the same for all the cells hitherto investigated. 

As an example, the dependence 1,(I.), in the blocking direction, for a cell coated with 
cadmium oxide is shown in Вела 


The relationships I, (V) and I, (V) for the same cell and for both the blocking and fore- 
ward directions of V are shown in Не. 3: | 


ano dark current 


] П | э — — —— 
10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 di 


[44] 


Fig. 2. СаО coated cell. Relation between photocurrent Jy and cell-current Г. 


300 200 100 0 100 200 300 


Fig. 3. CdO coated cell. Photocurrent Ir and cell-current J, against voltage V. 
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Fig. 4. SnO, coated cell. Photocurrent I; and cell-current 7, against voltage Ё, 


Fig. 4. represents the same dependences as in Fig. 3 for a specimen coated with 5лО , 

The following conclusions can be drawn from the graphs of Figs. 2, 3 and 4: 

a) “Dry” luminescent cells show a pronounced property of rectification. | 

b) For 1, considerably exceeding the dark current the dependence I.) is nearly linear: 
consequently L = cl, as in the electrolyte at constant oxide thickness. 

c) The shape of the curves in Figs. 3 and 4 suggests an exponential form of the fune- 
tions /,(V) and I,(V). 

In Fig. 5, log (I, +a) with suitably chosen values of the parameter a, is plotted against И 
for either cell. i 

As will be seen from Fig. 5, the relationship between log (I. +a) and V may be regarded 
as linear; hence, by eqs. (2) and (3), the dependence of the luminescence L on V can Ье 
represented by a formula of the type 


L = a(exp bV — 1). 


This is the relationship found for the electrolytic cells at constant oxide thickness. 


T 


log (Jet a) 
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Fig. 5. Log (I, + a) against V. I, in microamps. (1), СаО coated cell; (2), SnO, coated cell. 


E The L(V) 


and (И) curves from our experiments present the same general shape as 


hose obtained by D. A. Cusano and F. E. Williams (1956) in their investigations on the 


‚hotoelectroluminescence of ZnS:Mn, Cl films (Cusano 1955). 
For comparison, Fig. 3 from the paper of Cusano and Williams is reproduced here 


s Fig. 6. 


Both types of cells, viz., the АО» and 


the activated ZnS cells, exhibit the properties 


f rectification and of light emission in electric fields applied in the blocking direction. 
The only difference consists in the fact that, in the case of ZnS, the blocking direction 
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Fig. 6. Repreduction of Fig. 3 from the paper of Cusano and Williams (1956). 
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is from the semiconducting electrode, (TiO, 4-), to the metal, whereas in our cells the meta 
is on a positive potential with respect to the semiconductor. 

The curves of Figs. 3 and 4, in their general trend, also resemble the analogous curve: 
obtained by A. Fischer (1953) in his investigations of the electroluminescence of Zn( 
sheets, although the arrangement in Fischer's experiments was entirely different. 

Moreover our I,(V) curves are fairly similar to the analogous curve given by W. А 
Thornton (1959) in his paper on the electroluminescence of thin activated ZnS films. 


4, Concluding Remarks 


The fact of light emission from “dry’’ aluminium oxide films under the influence o 
strong electric fields, and the close similarity between our ЦИ) and L(V) characteristics and 
the analogous ones obtained in typical experiments on electroluminescence make it plaus 
ible to interpret the phenomena described in this article as electroluminescence effects 

On the other hand, the system Al/Al,O,/ semiconductor behaves in D.C. electric field: 
much in the same way as the system Al/Al,O,/ electrolyte at constant oxide thickness. 

Hence it would seem that, in either system, we are dealing with a phenomenon of true 


electroluminescence. 

It must be emphasized, however, that all our experiments were carried out in atmos: 
pheric air, which necessarily contains water vapour. 

This last remark is of importance, as some of our observations revealed a significan 
róle played by the humidity in the phenomena here described. 

The influence of the humidity of the medium will be the subject of a separate publica 
tion, 
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A continuity equation describing the motion of localized excitons in an anthracene crystal 
is proposed. The intermolecular transition probabilities are computed in dipol approximation 
with a semi-empirical adjustment as used by Dexter. Simpson’s experiment is discussed and 


his phenomenological theory is reexamined. 


Introduction 


A molecule in an anthracene crystal when excited may return to its ground state. 
leither by emitting a light quantum, or by transmitting the state of excitation to another 
olecule. A conversion of the excitation energy into lattice vibrations is also possible. 
his seldom happens, however, as the quantum efficiency of anthracene fluorescence 
is about 90% (Bowen, Mikiewicz and Smith 1949). A state of excitation travelling from 
nolecule to molecule is called an exciton. 

The transfer of excitons in anthracene was studied experimentally by Simpson (1956). 
In his experiment a thin layer of anthracene was deposited on a glass plate and covered by 
a layer of anthracene doped with naphthacene. This second layer worked as a detector of 
excitons because all the excitons entering it were trapped by naphthacene molecules and 
emitted as a green fluorescence, different from the fluorescence of anthracene. The double 
layer was illuminated by light falling through the glass plate, perpendicular to the layers. 
Let us denote the numbers of quanta per second: transmitted through the-anthracene layer 
in the absence of the detector, transmitted by both layers, and emitted by the detector, res- 
where J is the incident light intensity. If there were no exciton 
y of anthracene would be independent of the presence of the 


+ =a, would hold. Actually it has been found that «+a 
anthracene travel into the 


pectively by: al, Kl, г, 
transfer, the transmittivit 
detector and the equality a 
always exceeds x. It means that some of the quanta absorbed by 
detector and are emitted by naphthacene. The difference S=(a,+0,—@)I is the rate of 
transfer of excitons into the detector. The dependence of 5 on the thickness of the anthracene 
layer was determined by measuring & +a, —a, for 18 different specimens. 

(313) 
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To explain the results, Simpson assumed that the excitons are localized and that th 
flux of excitons, perpendicular to the surface, is given by the formula 


dn 


„en 1 
ردک‎ ( 


where D is a constant, n is the concentration of excitons, and the axis x is perpendicula 
to the layers. This formula implies that n fulfils a one-dimensional diffusion equation which 
for the stationary state takes the form i 

n d?n 

— — D— = е“ 

T dx? 3 ( 
where т is the mean life-time of an exciton and eis the absorption coefficient of anthracen: 
for the exciting radiation. For an anthracene layer of thickness d the following boundary 
conditions were assumed: 


dn 
DZ, for х=0, (3 
п. << () for x = d. (3b 


Equations (3a) and (3b) mean that the flux of excitons from glass into anthracene vanishe: 
and that the detector traps all the excitons which fall into it. This theory explains well tha 
experimental results when D is put equal to (460 А)?/т. 

In this paper a theory of exciton transfer is developed starting from the assumption thar 
each exciton is localized in one molecule at a time. This assumption cannot be strictly valid! 
for it is known that the Davydov splitting, observed also in the anthracene spectrum (Craig 
and Hobbins 1955), is caused by the delocalization of excitons; however, such a model car! 
be fairly easily handled and leads to some conclusions probably valid for real crystals 100) 


Theory 


Let us assume that the deposited anthracene layers are oriented with the cleavage 
plane (001) parallel to the surface of the glass. This result was experimentally obtained by 
Obereimov and Prichotko (1936), and confirmed by Craig and Hobbins (1955). The structure 
of anthracene crystals was studied by X-ray diffraction by Mathieson, Robertson and Sine- 
lair (1950). According to these authors the anthracene molecules in the crystal are arranged 
in layers, parallel to the (001) plane. From now on, the word layer is reserved for unimole- 
cular layers, and the whole anthracene layer is referred to as the plate. The distance between 
successive layers is taken as unit of length. 

Let us ascribe labels 1, 2, ... to successive layers, so that layer 1 touches the glass plate, 
and a layer labelled Ё is at a distance of k— 1 units of length from the glass surface. Let the 
index d label the first layer of the detector, then the layers 1, 2, ... ,d—1 consist of pure 
anthracene, and the layers d, d 4-1, ... belong to the detector. The number of excitons in 
the i-th layer will be denoted by n(i). In the stationary state the number of excitons pro- 
duced in a chosen layer together with those which reach it from other layers must equal the 
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umber of light quanta emitted from this layer together with the excitons transmitted from 


t to other layers. Consequently, assuming that the detector is very thick, we can write the 
Icontinuity equation as follows: 


А со д-1 
m + у n(i) pj — > n(k) ры = gle-*G-D, (4) 
k=1 k=1 Uic 2 и | 


where p; denotes the probability per unit time that an exciton present in the i-th layer 
ill pass to the k-th layer. 
In our notation the rate of transfer of excitons to the detector is given by the formula 


d—1 oo 
S(d) = 2 2 n(i) Ра: (5) 


"To obtain numerical values of the flux 5 which can be compared with experiment, we need the 
coefficients pj, and the numbers n(i). If the p;, are known, n(7) may, in principle, be found 
by solving equation (4). Thus our next task is to compute the transition probabilities р. 


Calculation of the transition probabilities pg, 


The probability of transition pj, is a sum of the probabilities of transition from a given 
molecule in the layer i — let us label it ij — to all the molecules kl from the layer Е. Thus 
the problem splits into two parts, first pj; ,;, the probability of transition between single 
molecules, is to be computed, next the result must be summed over all the molecules from 
the layer k (summation over J). It follows from symmetry that the result is independent of 
j and gives the required coefficient p;. 

In order to compute the transition probability for two molecules, let us assume that it 
is proportional to the probability computed by the perturbation method in the Born-Oppen- 
heimer approximation. The proportionality coefficient will be assumed to be independent 
of the relative orientation and distance of the two molecules (Golebiewski and Witkowski 
1959). Thus in the dipole approximation 

nis 

pij, = а ыш (6) 
‘where а is a constant and V;; j/R? is the interaction energy of the transition dipoles of 
molecules ij and kl, divided by the square of the dipole moment measured in units À x elec- 
tronic charge, thus |V;; „|?/ R9 is measured in À-$. The first optical transition of anthracene, 
which was active in Simson’s experiment, is short axis polarized (Craig and Hobbins 1955), 
consequently the transition dipole of the molecule points along its short axis. The magni- 
tude of the dipole moment influences only the coefficient a. 

Formula (6) gives the probability of transition to a single molecule. The summation 
over all molecules of the layer k which is now to be performed would be very laborious, 
and was therefore replaced by an integration. The result is 


Vi, 2 том я 
Уы _6x10 64-6, (7) 
1 
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The introduction of an integration instead of a summation increases the value of the sum. 
To get an estimate of the error thus introduced the part of the sum (7) corresponding to 
я | ; . . e . pe. 
molecules АГ parallel to the molecule ij was computed by direct summation for |k—il=1. 
The result is 2.5 x 10-6А-6 to be compared with the value 2.7 x 10-9 Д-6 obtained by inte- 
gration, thus the error is about 10%. Another source of error (Golebiewski and Witkowel 
1959) is the use of dipole approximation. A numerical estimate has shown that this error 
is of the same order as and of the opposite sign to the first one, thus both errors partly cancel. . 

From formulae (6) and (7), we get 


6 x 1079a b 


Pik secl - E sec, (8) | 
[6 Ta] 


The constant a is a semi-empirical parameter. For large intermolecular distances it 
1, екеш я Г 
may be estimated (Golebiewski and Witkowski 1959) from Dexter’s (1953) theory. In par- 
ticular if the average life time of the excitons т is known. a is given by the formula 
Ohl (9) 
а= ==. 
8n® E91? 
Here n is the refraction coefficient of anthracene, E is the transition energy from the region 
where the absorbtion and emission spectra overlap, and 5 is defined as follows: assuming 
J(E) and F(E) to be proportional to the numbers of quanta of energy Е, respectively emitted 


or absorbed by the molecule, and f (Е) Е = J F(E)dE = 1, we put 
= JKE F(E)dE. (10) 


The parameter 5 causes the most difficulty when applying this theory. It depends 
critically on the shape of the spectra, and consequently on the method of preparing the 
crystal (Birks and Cameron 1959). What is more, it is not at all clear whether p should 
depend on 5 for aromatic hydrocarbons. Northorp and Simpson (1956) deny this on the 
basis of experiments with the quenching efficiency. It should be borne in mind therefore 
that the coefficient a, computed from (9). is at best an estimate. 

For numerical computations the values: т = 6.4 x 10 sec (Birks and Cameron 1949), 
n = 1.6, E = 3.1 eV (Bowen 1949), and S =2 x 10-2 eV-1 were assumed. The value of n is 
estimated from the values for related compounds (Handbook). The value of S is a crude 
estimate from the graphs given by Bowen (1949 — with no scale of intensities), and Birks 
and Cameron (1959 — without the absorption spectrum). Substituting these values 
a = 4.4 x 1015 Ав sec-1 was obtained, which gives 


b = 2.7 x 1010 вест, (11) 


There is one more problem to be discussed in connection with formula (8). Up to 
now shielding has been completely neglected. If we assume on the contrary ideal shield- 
ing, only p, ,,, can be different from zero. Actually it should be less than computed by 
formula (6), but not much less, as the sum (6) converges rather quickly (the largest term 
yields by itself about 10% of the whole sum). In the following the decrease of py p+ Caused 


by shielding is neglected. In the usual theory of sensibilized luminescence shielding is not 
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cluded, but it will be shown in the following section that this assumption is necessary 
o derive the diffusion equation from equation (4). For the sake of brevity the models with 
nd without shielding will be referred to as models (a) and (b) respectively. Both models will 
ye dealt with in the next section. | 


Discussion of the continuity equation 


The system of equations (4) may be solved easily if model (a) is adopted. In this case 
he equations can be written in the form 
n(i) Е дыд J 
o Ръчна С 1) — 2т(0) + n vs Nele (11) 
t=1,2,...,.d—1 


here n(0) and n(d) are defined by the conditions 


раа) RO) =0 (0) = 0 (12) 
ssuring the equivalence between systems (11) and (4). 
It can be verified by direct substitution that the solution of this equation is 


etle? З 
n(i) = — pac cC om 


1— Api, rey T sinh? = 


(13) 
where the constants C, and C, should be computed from conditions (12) and 7 is the positive 
root of the equation 


Ар, ;,4t sinh? 4/2 = 1. (14) 
Actually в is about 7.5 x 107? and p; 1722250 (from experiment), thus the operators sinh 
and e* may be dropped and (13) passes into the solution of the diffusion equation found by 
Simpson with 


Dan (15) 


^j i31 
An alternative method of deriving this result is to assume equation (11) as valid for a conti- 
nuous variable x and to put in (11) and (12) differentials instead of differences. With the 
assignment (15) equation (2) with boundary conditions (3a) and (3b) is thus obtained. 

Formula (15) has been derived taking the distance between neighbouring layers c= 
—9.19 А (Mathieson, Robertson and Sinclair 1950) as the unit of length. Changing into 
angstroms, we get from (8) and (15) D —(121À)?/v to be compared with the experimental 
result D=(460A)2/r. Let us note further that for this model the flux of excitons falls expo- 
nentially with increasing thickness of the anthracene plate (Simpson 1956): 


g 1 SV Pr ened, (16) 


The solution of system (4) with the assumptions of model (b) is much more difficult. 
In this paper only an estimate, valid for thick layers, is given. This estimate, however, shows 
clearly the difference between the implications of models (a) and (b). 


318 


Whatever the distribution of excitons in anthracene the flux of excitons streaming into 
the detector is larger than the flux бу computed on the assumption that all excitons are 
localized in the molecular layer labelled 1. For 5, we get from formula (5) 


3 жы eue 1 | 
Nr е^ age 12) 


where N is Ше total number of excitons in anthracene. In unit time anthracene absorbs 
I(1—e~*4-) and emits N/r-- S quanta. Thus for the stationary state 


№ = h(1 — e-*4-») — Sz. (18) | 
Substituting into (17) and solving the inequality S, < 5 for S we get 


br 
z3;— (1 — ere@-D), 19): 
затоа р (19) 
We see that model (b) predicts a much slower decrease of the flux S with the thickness d of 
the anthracene plate than model (a). 


Conclusions 


The assumption that each exciton is localized in one molecule at a time implies equation 
(4), which is not equivalent in general to the diffusion equation proposed by Simpson. 
The subsequent caleulations depend on whether shielding is taken into account. 

In model (a) it is assumed that the shielding is perfect, in the sense that excitons can 
be transferred between neighbouring molecular layers only. In this case a diffusion equation 
for the concentration of excitons is obtained and an exponential decay of the exciton flux 
follows. The transmission coefficient computed by Dexter's theory is about an order of 
magnitude smaller than required to reproduce the experimental results. 

In model (b) shielding is neglected. In this case equation (4) does not reduce to a diffu- 
sion equation. The qualitative difference between models (a) and (b) is easily seen for thick 
anthracene layers. Here model (b) predicts a decrease of the flux, but only with the third 
power of the layer thickness. If this effect were really observed, it would shed a new light 
on the band theory of excitons where the applicability of the diffusion equation is assumed 
(Heller and Marcus 1951). 

The theory of this paper was developed neglecting the transfer of excitations by reabsorp- 
tion of fluorescence. Indeed, Simpson verified that in his experiment this effect is unimportant. 
For thick plates, however, it must be included, as proved by the observed dependence of 
the decay time of fluorescence on the size of the anthracene sample (Birks and Little 1953). 


REFERENCES 


Birks, J. B., Cameron, A. J. W., Proc. Roy. Soc., 249, 297 (1959). 
Birks, J. B., Little, W. A., Proc. Phys. Soc. [London], 66, 921 (1953). 
Bowen, E. J., The Chemical Aspects of Light Oxford 1949 ра162: 


319 


wen, E. J., Mikiewicz, E., Smith, F. W., Proc. Phys. Soc. [London], 62, 26 (1949). 
aig, D. P., Hobbins, P. C., J. Chem. Soc., 2309 (1955). 
exter, D. L., J. chem. Phys., 21, 836 (1953). 
olebiewski, A., Witkowski, A., Roczniki Chemii, 33, 1443 (1959). 
ndbook of Chemistry and Physics, Ohio Cleveland 1956 ed. 38. 
eller, W. R., Marcus, D. L., Phys. Rev., 84, 809 (1951). 
athieson, А. McL, Robertson J. M., Sinclair, V. C., Acta cryst., 3, 245 (1950). 
orthorp, D. C., Simpson, O., Proc. Roy. Soc., 234, 136 (1956). 
bereimov, I. W., Prichotko, A. F., Phys. Zs. Sow., 9, 48 (1936). 
impson, O., Proc. Roy. Soc., 238, 402 (1956). 


1. ХХ (1961) ACTA PHYSICA POLONICA ГА 


POSITRON RADIATION OF Tm’? 


By J. Zyticz, Z. PnEIBISZ 


Institute for Nuclear Research, Polish Academy of Sciences, Warsaw, 


S. CHOJNACKI, J. WorowskKki 


Institute of Experimental Physics, Warsaw University, Warsaw, 


Yu. NORSEEV 


Joint Institute for Nuclear Research, Dubna, USSR. 
(Received October 20, 1960) 


The Tm!®% positron spectrum has been examined with a long-lens beta spectrometer fitted 
with helical baffles for separating the electrons from the positrons. Two components were found 
of maximum energy 1920--20 keV and 1219--40 keV and relative intensities 1: (0.092-1- 0.018). 
'The total number of positrons for a decay amounts to about 1%. The decay scheme of Tm!** 
and the applicability of the theory of non-axial nuclei is under discussion. 


I. Introduction 


| Data on positron radiation of Тт! (Ty). = 7.7h) can be found in papers by Wilkinson 
nd Hicks (1949), Boskma and De Waard (1959), and Wilson and Pool (1960). According 
to the first of these, the maximum energy of the positron spectrum amounts to 2.1 MeV 
(about 0.595 positrons per decay). Boskma and De Waard claim to have found two compo- 
nents of energies 450 + 80 keV and 1.300 100 keV (relatively 1.8% and 1.3% positrons 
per decay). Lastly, Wilson and Pool, in a short note announcing the publication of a more 
detailed paper, report that according to their measurements the maximum energy amounts 
to 2090 + 40 keV. In each of the three works pure! Tm!99 was produced and then in- 
vestigated. 

The aim of the present work was: 1. to measure precisely the energies and intensities 
of the components of Tm1$6 В+ radiation and thus to solve the inconsistency observed in 
the data so far published, 2. to discuss the decay scheme, 3. to consider the possibility of 


applying the theory of non-axial nuclei. 


1 Not in equilibrium with Yb. 
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II. Source. Apparatus. Measurements 


Isotope Tm!955 was produced in the spallation process of tantalum nuclei o b» 
660 MeV protons from the ЛУК synchrocyclotron. The rare earth elements were chemically 
separated from the tantalum target irradiated for four hours. Next, the rare earths were 
separated according to their atomic number by a chromatographic method (Preobrazhensks 
et al. 1956). 

Tm166 is produced directly in a spallation reaction or by decay of Yb1® (Т, = 57 h)} 
so that either the thulium fraction or ytterbium fraction may be used for measurements: 
It appears that the use of the ytterbium fraction is more advantageous, as a greater activity 


Тт 


Jntensity 
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Fig. 1. Results of additional separation of the ytterbium fraction. 


slower to expire is obtained. Yb166 does not e 
does not hinder the observation of Tm166, 


isotope present in this fraction: Yp169. 


mit positrons and so its presence in the source? 
The same can be said of the second long-lived| 


To avoid any possible contamination, the 


ytterbium fraction was subjected once again! 
to the chromatographic process shortly 


after its separation. The results are shown in Fig. 1 
where the activity of drops leaving the column is plotted as a function of the drop number.. 
Three peaks are visible: the weak peak of lutecium — the contamination, the basic peak: 
of ytterbium, and the thulium peak — the Tm!6 produced from Yb166 during the time: 


interval between the first and. second chemical operation. In the preparation of the source, . 


99 
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ops of the ytterbium peak were used. The drops were evaporated on a thin aluminium foil 
136 mg/cm?) moistened with an insulin solution. When the measurements were started? 
‚hours after chemical operation, Tm!6 was already almost in equilibrium with Yb166 
id its activity amounted to about 0.5 mC. 
A spectrometer with a long magnetic lens and baffles for separating the positrons from 
е electrons (Chojnacki et al. 1960) was used for carrying out the measurements. The mica 
indow of the G.M. counter was 1.3 mg/cm? thick and did not cause any distortion of the 
Jectrum for particles with energy exceeding 170 keV. 


150 


мее) 


0 500 1000 1500 2000 
E(keV) 


Fig. 2. Energy distribution of Tm!$9 positrons. 


Three measurements of the positron spectral distribution were made. Each measure- 
ent took about four hours, so that small corrections for the expiration of activity with time 
ere necessary. Аз an example, Fig. 2 shows the results of one measurement. 

It was verified during a period of thirteen days that the activity of the positron radiation 
minished throughout the entire spectrum in equal proportion with a half-life of 57h. 
his result was in agreement with the assumption that the radiation under observation 
is due to Tm196, the latter being in equilibrium with Yb". 

In addition to the positrons, internal conversion electrons from the Tm!96 were also 
easured. Part of the electron spectrum is shown in Fig. 3. 


III. Results 


Fig. 4 represents the positron spectrum in the form of a Kurie plet. With the assumption 
an allowed shape of spectrum the analysis yielded two components with maximum пеш 
1.9 MeV and 1.2 MeV and of relative intensities of 10:1. The presence of a third, less 
tense component of maximum energy lower than 1 MeV cannot be ruled out. 


2 After moving the isotopes to laboratory in Warsaw. 


324 


16000 к- 459.3 


- 


14 000 f 


12000 


T 


10000 + 


Counts per 2 min 


13 14 15 16 17 
Coil current in amps 


Fig. 3. Part of the internal conversion electron spectrum of Tm!*5, 
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Fig. 4. Kurie plot of the Tm!** positron spectrum. 


The computation of log ft for both components required the determination of tl 
mean number of positrons per Tm16 decay, This was deduced from data on the conversi 
radiation of this isotope. 


l. The ratio of intensities of the Pj component (Emax = 1.9 MeV) and of the K-co 
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reion line of the 459.3 keV ? transition was measured: 
ВЕ: ері ао з = 1: (0.079 + 0.0095) 


e K-459.3 line was chosen for the measurement because it is well resolved from the 
ner conversion ‘lines and its activity is not excessive for measurement with a G.M. 
unter. 

` 2. On the basis of a paper by Gromov et al. (1959), the absolute intensity of the K-459.3 
e was computed at (6.56 + 0.94) x 10-4 electrons per decay. 

` 8. From results 1 and 2, the absolute intensity of the Bi component is found to equal 
3 + 1.6) x 10-3 positrons per decay. 

` The ratio of the probabilities of electron capture and of positron emission (e/B*) was 
lculated using data assembled by Wapstra et al. (1959). These are exact for allowed transi- 
s and may serve as a first approximation for first forbidden transitions (probably the 
nd we dealt with in the present investigation). 

` The final results of all the measurements and computations are assembled in Ta- 


e І. 


TABLE I 
| я Я 
Relative Number of positrons 
Component E (keV) | des per decay | log ft e/B* (theor.) 
Bt 1920 + 20 1 (8.3 + 1.6) x 101 96 7.9 ~ 4 
Bt 1219 + 40 0.092 + 0.018 (7.6 + 2.1) x 10? % 8.1 ~18 


IV. Control measurements 


Element fractions of the rare earths are separated by the chromatographic method 
ithout using a carrier; hence, sources prepared from these fractions are of very small 
tickness. Nevertheless, the spectrum in the low energy region can still be distorted due 
electron scattering. In order to estimate the amount of distortion due to apparatus factors 
| our case, a source was prepared of a neodymium fraction from the tantalum target and 
je positron spectral distribution of Pr49 obtained from №10 decay was measured. In 
-urie coordinates a straight line was obtained beginning from maximum energy (2362 keV) 
own to an energy of about 350 keV. This led to the following conclusions: 1. the Pr!# 
ositron spectrum above 350 keV is simple and 2. above this energy electron. scattering 
oes not cause distortion of the spectral distribution. Conversely, in the case of thulium, 
nall deviations from the straight line in the Kurie plot of the f; component (Fig. 4) may 
lready be observed for energies of about 600 keV. 

The positron spectral distribution of pure Tm199 separated from an ytterbium fraction 
as also investigated. The analysis of the Kurie plot confirmed the existence of components 
+ and fj. No traceswerefound of a third component; this, however, may be due to low 


tatistical accuracy. 


3 Here and subsequently the enegies of Tm!#® y transitions are quoted according to Jacob et al. (1960). 
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V. Discussion of the decay scheme 


The difference between the maximum energies of the components f, and f, amoun 
to 701 +45 keV. An analysis of the Er!#® excited levels scheme (Fig. 5) from the point « 
view of possible ways of positron decay leads to the conclusion that two pairs of levels a: 
distinguished, namely 80.6 (2*) keV and 787.1 (27) keV (the difference in excitation energ 
is 706.5 keV), and 265.1 (47) keV and 957.2 (47) keV (with a difference in excitation епего 
amounting to 692.1 keV). 


Tm 166 (778) 


+ а 7 
9572 4 — ера иа: 
ER «а : / > ® 
860.6 3* dA / Pr x 
787.1 : X 
y 
| | 
545.3 | 
265.1 
80.6 
0 Жо че PS ET 
keV 


Fig. 5. A possible decay scheme of Tm!99, Scheme and energies of excited levels according to Jacob et al (1960) 


The assumption that the components Û, and P» excite 80.6 keV and 787.1 keV levels 
is in agreement with what was known on isotope Tm166 previous to Wilson and Pool: 
publication (1960), namely 1. that the theoretical aspect (comparison with Tm1® and Dy163 

. € с Ь . r 1 
suggests spin 2 or 3 and odd parity for the Tm!99 ground state, and 2. that the experimental 
data distinguish the first of these values (Gromov and Dzhelepov 1959, p. 65—66). Indeed 
if the Tm! ground state is (2°), then it seems obvious that positron transitions will lead 
primarily to levels (2*). Thus we have the first forbidden transition. This is also suggested 
by the log ft value (see Table 1). The decay scheme of this variation is shown in Fig. 5. 
The transition energy 0-3028 keV does not appreciably differ from the 2700 keV value 
obtained from the semi-empirical equation by Cameron (1957). Levy’s (1957) equation yields 
the distinctly smaller value of 2230 keV. 

| Wilson and Pool (1960) report that according to their analysis the principal positron 
radiation component excites the level 265 (4*) keV. In this case the transition energy Q would 
amount to 3207 keV. 
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VI. Possibility of applying the theory of non-axial nuclei 


According to Davydov (1959a), the nucleus of Ег188 is one of those not possessing the 
dal symmetry. If this is the case, the first three excited levels (normal) as well as the next 
ree (abnormal) have a rotational character (Fig. 5). Rotational nuclear excitation has no 
preciable effect on the internal state of the nucleus. In other words, the wave function 
escribing the internal state of the nucleus is the same for all rotational levels. This allows 
s to generalize the results obtained by Alaga et al. (1955) and to write the equations for 
e relative intensities of the f transitions in the case when one transition leads to a “normal” 


vel and the other to an “арпогта]” one (Davydov 1959b). These equations may in principle 
e applied in the case of Tm1® positron radiation. Nevertheless, some difficulties subsist. 

1. It is not certain what the spin of Tm?® is in its ground state. 

2. In general, it should be assumed that the ground state of Tm!99 is a superposition 
f states corresponding to different values of the quantum number К. However, there is no 
heory for odd-odd nuclei (such as Tm!99) which would permit us to evaluate the contribu- 
Лоп of the particular states. 

3. Davydov's (1959b) equations are applicable to transitions of a fixed multipolarity L, 
hereas in the case of forbidden transitions of the first order L=0, 1 and 2 should be exam- 
ned. It is not possible, however, to predict the ratios of the corresponding matrix elements, i.e. 
o predict the contributions from the different multipolarities. 

Taking into consideration the foregoing difficulties: 
| 1. It is assumed that the Tm! ground state has spin 2 (Ше most probable value) and 
odd parity. 
| 2. It is accepted that the ground state of this isotope is a pure К —2 state (probably 
the contribution of other states is much lower). 

3. Possible first forbidden transitions are taken into consideration: f' — to a lower 
level (2^), B" — to a higher level (27), and p" — to level (3+) (Fig. 5) and the different 
cases of L—0, 1 and 2 are examined. 

The y parameter, determining the deviation of the nucleus from axial symmetry amounts 
to 12.7 for Er!6*, Knowledge of this number enables us to give the wave functions of the 
three excited states already mentioned (Davydov and Filippov 1958). These are the following: 


Фл = 1.000 Dy + 0.017Ф,, 


Poom = 0.017 Da + 1.0009, 
= Da 


W31m 


Фк is the wave function of a state with spin J and quantum number K. 
The ratios of the reduced probabilities of transitions f and f" to the reduced proba- 
bility of transition В” are given in Table II for different L.*) 


4 [n the case of L —2 equal intrinsic . matrix elements of K =2—K 22 and K —2— К = 0 


transitions were assumed to show, that there is no K-forbidness. However the real intrinsic matrix 


elements may not be equal. 
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TABLE II 


ratio of reduced probabilities | L=0 ESI JL 


B’ (L; 2— 21) 


nt рар Ses 2.9x 1074 2.9x 1074 1.0 
B" (Г; 2> 22) 


B” (L; 2—31) 
B” (L; 2-> 29) 


| 0 0.50 | 1.8 


The ratio of the reduced probabilities (Bi|Bs),,, found experimentally for Ше compo- 
nents f, and f, amounts to 1.44. If it is assumed that the f component experimentally 
observed is a superposition of f" and f" transitions and the f, and 8" components are 


identical, then 
m >|-- << 1.44 
Е | ехр В, | exp 


Thus, with the all foregoing assumptions it is not possible to obtain a sufficient degree of 
consistency between the theoretical and experimental data. 


VII. Supplementary note 


The present paper was presented at the Conference on Neutron Deficient Isotopes at 
Dubna in June 1960. During this conference it became known that Tm199 positron radiation 
had also been analysed by Brovtsyn, Gromoy, Dzhelepov, Dmitriev, and Nikitin. Their 
results (as far as we know of a preliminary character) were presented at the Conference on 
Nuclear Spectroscopy in Moscow in January 1960. Brovtsyn and co-workers found two p* 
components: 2020 +20 keV and 1150 +100 keV of relative intensities 1:0.17 and stressed 
the probability of the existence of a third component: 220 keV. 

The authors wish to express their thanks to Docent Z. Wilhelmi, Head of the Nuclear 
Physies Laboratory of the Institute for Nuclear Research in, Warsaw, for his kind interest. 
in the present investigation and to Miss I. Jarosiewicz and Mr Z. Brynikowski for their 
help in the measurements and computations, The authors feel especially indebted to Pro- 
fessor B. Mottelson for his invaluable comments on the discussion. 
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ОМЕ EXAMPLES OF INTERACTIONS OF VERY HIGH ENERGY 
ROTONS WITH HEAVY NUCLEI OF PHOTOGRAPHIC EMULSIONS 


By J. BARTKE, Z. CzacHowska, R. HoryNskı AND K. RYBICKI 


Cosmic Ray Department of the Institute of Nuclear Research Cracow — Poland 


( Received October 24, 1960) 


Three high energy jets (Ep,jm> 10" eV) produced in collisions with nuclei of photographic 
emulsion are described. In spite of the fact that they are probably produced in central collisions 
of nucleons with heavy nuclei they show double maximum angular distributions in contradiction 
with hydrodynamical model. 


Introduction 


In one of the previous papers of our laboratory (Bartke et al. 1960) it has been pointed 
out that there exists a double maximum angular distribution in events for which we have 
evidence, that they are central collisions of nucleons with heavy nuclei. Since this problem 
seems to be important from the point of view of different models of multiple meson produc- 
tion in very high energy collisions we would like to describe three events of this type. 


Experimental 


In the scanning for nuclear interactions of cosmic particles three jets of comparatively 
high energy with great numbers of both evaporation and shower particles have been found: 
jet I: 26 +47p, jet П: 18 +-41р, and jet Ш: 15 +78p 1, The events are characterized by a great 
multiplicity and high excitation of the target nucleus and can be interpreted as central 
collisions of nucleons with heavy nuclei (Ag or Br). 

Jets I and II have been found during a scanning for jets via electromagnetic cascades 
in a stack of Ilford G5 emulsions irradiated in the Po valley in 1957, jet Ш was found in 
area scanning in a NIKFI-R emulsion stack irradiated near Moscow in 1958. 

In all these events we have measured angular distributions of minimum, grey and 
black tracks. The results of these measurements and their analysis are presented in sec. 2 


and 3. 


1 These jets have been referred to in the paper of Gierula et al. (1960), as 171K, 168K, and 200K res- 


pectively. 
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ТАВГЕ 1 
Laboratory angles of secondary particles with respect to the jet axis for jets: I, II, III. 
Jet I 
ps log tan 0 Kar log tan 0 er log tan 0 
1 5.825 17 3.189 33 2.987 
2 4.404 18 3.190 34 1.056 
3 4.405 19 3.233 35 1174 
4 4.586 20 3.236 36 1.321 
5 4.691 21 3.293 31 1.355 
6 4.781 22 3.309 38 1433 ^ 
7 4.826 23 3.457 39 1.559 
8 4.859 24 3.770 40 1.573 
9 4.871 25 2.242 41 1.712 
10 4.949 26 | 2.286 42 1.746 
11 4.972 27 2.357 43 1.764 
12 3.123 28 2.361 44 1.860 
13 3.125 29 2.457 45 0.061 
\ 14 3.125 30 2.648 46 0.188 
15 3.151 31 2.879 AT 0.247 
16 3.185 32 | 2.962 
Jet II) 
pes log tan 0 s log tan 0 e log tan 6 
1 4.532 15 3.447 29 1.099 
2 4.634 16 3.462 30 1.107 
3 4.681 17 3.613 31 1.110 
4 4.699 18 2.204 32 1.202 
5 4.748 19 2.249 33 1.231 
6 4.892 20 2.263 34 1.259 
7 4.964 21 2.290 35 1.317 
8 3.041 22 2.536 36 1.339 
9 3.079 23 2.539 37 1.723 
10 3.204 24 2.543 38 1.819 
11 3.322 25 2.637 39 1.822 
12 3.362 26 2.789 40 1.983 
13 3.415 27 2.885 A 0.244 
14 3.447 28 2.952 


*) Four tracks not presented in the table (0,—60,—6.9 x 10-4, 05=6.2 х 10? and 0,—6.9 x 10~rad.) 
are considered as electron tracks. A definite determination of the nature of the particles is not possible because 
of great track density in the inner cone and the development of the electromagnetic cascades in further plates. 
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TABLE I (continued) 


Jet III 
Eb log tan 6 Bi log pam 0 лы (pt y" 

1 3.342 26 2.681 51 1.134 
2 3.362 21 2.690 52 1.134 
3 3.602 28 2.699 53 1.146 
4 3.699 29 2.128 54 1.218 
5 3.708 30 2.785 55 1.281 
6 3.740 31 2.803 56 1.295 
1 3.763 32 2.813 57 1.305 
| 8 3.813 33 2.823 58 1.308 
9 3.851 34 2.826 59 1.337 
00! 3.860 ` 35 2.833 60 1.337 
п 3.886 36 2.900 61 1.340 
12 3.898 37 2.914 62 1.342 
13 3.078 38 2.919 63 1.378 
14. 3.996 39 2.935 64 1.515 
15 2.033 40 2.944 65 1.593 
16 2.041 41 2.954 66 1.613 
17 2.051 42 2.969 61 1.671 
18 2.090 43 2.978 68 1.681 
19 2.111 44 1.009 69 1.683 
20 2.290 45 1.025 70 1.686 
21 2.371 46 1.029 za) 1.825 
22 2.484 47 1.045 2 1.868 
23 2.491 48 1.057 73 1.939 
24 2.648 49 1.061 74 0.455 
25 2.658 50 1.097 75, 0.735 


Three tracks are not included since they are emitted backwards. 


We have also scanned for secondary interactions. In jet I the scanning has been per- 
formed only in the narrow cone (volume containing 23 tracks with the total length of ~110 cm) 
and no secondary interactions have been found. In jet П ме have scanned an area containing 
32 tracks with the total length of —260 cm. Three secondary jets was found i.e. 4 +-9р, 
with y, —13, produced in the distance of 49 mm from the primary interactions by a singly 
charged secondary emitted with an angle 0, —14'; (transverse momentum ру —1.4 BeV/c), 
and also jets 14 +(28 +2)p, 2.5<y,<3.1 and 26 +(>20)p, 10<y,<20. A better determina- 
tion of the track number and the y, — value was not possible because of poor visibility 
and great distortion in this part of the stack. In the jet III a secondary interaction was found 
already in the distance of 196 from the primary interaction point. This secondary interac- 
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tion 12+6p or n was initiated by one of the tracks of the narrow cone. Target diagrams in 
the distance 600u, 10004 and 1400и from the primary interaction enable to distinguish 
between the tracks from the secondary interaction and those from the primary. 


2. Angular distribution of secondary particles and the determination of the primary energy. 


For each jet we made a target diagram by measuring the coordinate y and z with respect | 
to one selected track. The polar angles of the tracks of the narrow cone with the jet axis: 
were determinated by measuring the distances of the points on the target diagram from the: 
centre of gravity of this points. In the diffuse cone these angles were determined from 
the projected angle g (on the plane of emulsion) and the dip angle 9 using the formula: | 
sem 0; = sem p+ cos Ü cos Jy + sem (0 — 00), 9, — being the dip angle of the primary. . 
The results are presented in table I. 

From the angular distribution the Lorentz-factor of a system in which there is for- 


ward-backward symmetry can be determined as y, = log ctg 0;, 0; — being the polar: 
angle of the i-th particle with respect to primary particle direction (jet axis). 

For nucleon-nucleon collision we obtain the primary energy in the laboratory system 
from: E, = 2у2 Ме? (y, > 1). 

In collisions of nucleons with nuclei we must assume a type of mechanism in order 
to evaluate the primary energy. If we assume that the primary nucleon interacts only with 
a group of the nucleons of the target nucleus (tunnel model) we obtain for the energy: 
E,=2ny?Mc? where п is the number of nucleons in the tunnel. In our events using this 
model (collisions with Ag or Br) we can assume average n=4. 

If the angular distribution is presented in terms of dN/d (log tan бү) vs. x=log tan 0, 
we can use as the measure of the anisotropy the dispersion 


ME y^ (log tan 0; — log tan б) 
п — 1 


For the isotropic distribution this dispersion should be о—=0.39. Anisotropic distributions 
are characterized by greater o-values. 


7 


In the table II the E,, y, and о -values for the events described are presented. 


| TABLE II 
jet | type | 7 | E, = 2у2Мез | E,—2nyiMe | с 
| | 
I 26+47p 83.1 1.3 x 1013 eV 5.2x 1013 eV 1.25 
I 18+41р 58.1 6.4х 1012 eV 2.5 x 1013 eV 1.10 
IH 15+78р 15.3 4.4x10%eV |  L8x109eV 0.71 


On Figs. la, b, с, and 2a, b, c, we present differential and integral angular distributions 


of the jets in the coordinates x — log tan 0, and y = log I where F(0) is the fraction 
of number of tracks with angles smaller than 0. 
It can be seen that the differential angular distribution shows strong anisotropy and 


а double maximum structure which corresponds to a plateau in the integral dostribution 
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Fig. la, b, c Differential angular distributions for jets a — I, b — II, c— Ш. The continuous curve represents 


he best fit of Gussian curve to the experimental data with o = 1.25; o = 1.09 and o= 0.71 respectively. 


336 


(og 
1 
-5 -3 
log tan8 T 
= = 
. rat 
a 
; D: 


EMEN SI ШЫ 


log tan 8 oo! 


Fig. 2a, b, с Integral distributions for jets a— I, b — II and c— III in Duller-W. 


alker's coordinates. е 
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3. Angular disiributions of black and grey tracks 


We have measured the angular distributions at black and grey tracks for all events. 
ey have not significant deviations from the isotropic distribution, what may be seen in 
в. За, b, c. The estimated total energies of black and grey tracks are: jet 1: 900 MeV 

II: 1400 MeV and jet III: 500 MeV. | | 


0 Ju "il 0 *1 


с. За, b, c Differential angular distributions of the black and grey tracks for jets а — I, 5 — П and c— III 
coordinates cos бу. The continuous and dotted lines correspond to black and- grey tracks respectively. 


Comparison of the angular distributions of the secondary particles with the predictions 


of the hydrodynamical model 
The hydrodynamical model (Belenkii-Landau 1956, Milekhin 1958) has predicted 
е normal shape of angular distributions with the maximum for CM-angle 0 — 5, both 


r nucleon-nucleon and nucleon-nucleus collisions. 

In the original paper in which the double-maximum distributions had been described 
"iok et al. 1958a, 1958b) only nucleon-nucleon collisions were considered. The two-centre 
odel has been introduced for explaining this shape of distributions for n-n collisions. 
aking into account the possibility for the existence of this type of angular distributions 
je opinion has been expressed. (Feinberg 1960) that roughly speaking there are two types 
f nucleon-nucleon collisions. There are collisions with complete transfer of energy (head-on 
‚llisions, inelasticity K=1) and collisions with incomplete transter of energy (peripheral 
Ilisions, inelasticity K<1). According to these ideas the hydrodynamical model in its 
resent form could be applied only to the first type of collisions. But especially favourable 


jnditions for applying of this model are fulfilled for central collisions of nucleons with 


eavy nuclei. 
According to hydrodynamical model these collisioı 


uclear matter of the tunnel which the nucleon shears 


1s correspond to collisions with the 
off the target nucleus. The number 
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of nucleons in the tunnel is characterized only by the length of the tunnel which is consider 
as without any structure. The differential angular distribution can, according to this theoı 


be well described by a Gaussian curve: dN/dx~ e-*?.: where L = 0.56 In E/M+ 1.6 


x In 4- 1.6. Thus L is a function of the tunnel length and the primary energy. 
1 


The analysis of the experimental facts does not confirm the predictions of this mod 
In the paper Bartke et al. (1960) a statistical analysis of nucleon-nucleus collisions w 
performed in a large energy interval and it was found that for jets with a primary ener; 
greater than 1012 eV and with a great number of evaporations tracks and shower particl 
there is a strong anisotropy and a double maximum structure. 

Taking into account the importance of the existence or not existence of the two-maxim 
-jets in central collisions with heavy nuclei, we present here in details the three cases fro 
our laboratory which may be considered as central collisions of nucleons with heavy nucl 
of photographic emulsion. From these cases I and IIT was not included in the paper « 
Bartke et al. (1960) and II was carefully remeasured. 


| 
one track 


- оо О * oo 
a b с 

Fig. 4a, b, c Differential angular distributions for jets a-I, b-II and c-III in equal area intervals of the Gaussid 

function. The Gaussian distribution is given by the dotted straight line. 


We have analysed the deviations of the distribution from the normal one in three ways 
1°, On Fig. 4a, 6, c, the angular distributions of the analysed jets are presented i} 
terms of histograms in coordinates in which equal intervals on the x 


-axis correspond tj 
equal areas of the normal distribution. The 


normal distribution curve is in this coordinate 
a straight line parallel to the x-axis. The deviations from the normal 


by means of the Pearson’s X^-est are presented in table III. 
29. In the paper of Gierula et al. (1960) 
deviation predicted by the two-centre mode 


distribution estimate: 


the authors introduced a new measure of th: 
l. If we divide the angular distribution in th: 
C.M. system into two parts by the value x = + 0.6740 the number of tracks outside thi: 
interval N, (i.e. those with |x| > 0.6740) and the number inside N; (|x| < 0,6740) should 


be equal for the normal distribution. In the сазе of the double maximum distribution pre 


aW N 
dicted by two-centre model the N, is greater than N,. Thus the value D = MN = N 0 
| №, + №; П. 


was proposed as a measure of this deviation, since it should be zero for the normal distri: 
bution. In the case of any deviation towards a symmetrical double maximum distribution 
this value is 0< D «1. The level of the significance of the deviation can be computee 
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: 4° = D?n,. The D values and the significance levels for the jets presented in this 
er are given in the 3° and 42 column of table III. 

39. The application of the Kolmogorov-Smirnov test (Kolmogorow 1953, Smirnow 
9) is a not quite correct procedure since the test in its present form cannot be applied 
comparing an experimental curve with a theoretical curve with parameters estimated 
m the experiment. The application of the tables available gives too low values for the 
iation. Nevertheless we have calculated these values according to this test (4% column 


le III). 


TABLE III 
Deviation computed by Deviation computed Deviation computed by 
means of the y?-test 2 from x?(D) a е опори 
rov-Smirnov test 
1 2 3 4 

1.5 st. dev. 0.53 3 st. dev. 2 st. dev. 

~1 st. dev. 0.41 2.5 st. dev. 1 st. dev. 

2.5 st. dev. 0.07 ^1 st. dev. 1 st. dev. 


5. Conclusions 


l. It has been confirmed on further events found in our laboratory that in central 
isons of nucleons with heavy nuclei for which there is a strong anisotropy we observe 
ouble maximum structure significantly different from the normal one. 

2. This fact is in disagreement with the hydrodynamical model. 

3. This fact would be in agreement with the hypothesis of the superposition of peripheral 
llisions each described by the two-centre model, if the "core" for "central" nucleon- 
ucleon collisions was much smaller than the nucleon radius. 

It is a great pleasure for us to thank Professors M. Miesowicz and J. Gierula for continuous 


idance, permanent discussions and hints during progress of this work. 
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ON THE SHAPES OF ALLOWED. BETA SPECTRA 


By Br. KucHowicz* 
(Joint Institute for Nuclear Research, Dubna, USSR) 


(Received October 26, 1960) 


The present paper brings a discussion of the small deviations observed in the beta spectra 
of 32P, 9°Y, 1141} and Na. No explanation thereof exists within the framework of the conven- 
tional theory, even assuming the usually dominating nuclear matrix element to be suppressed. 
Two alternative explanations of the phenomenon are discussed, both of which are tentative 
and of a somewhat surprising character. The one ascribes the effect to G-nonconserving interac- 
tions, whereas the other hypothesis attributes the deviations to a virtual intermediary chiral 


spinless boson. 


1. Small deviations in beta spectra 


It would seem that the foundations of beta decay theory are rather well established, as 
e Universal Fermi Interaction explains a considerable number of the phenomena. There 
main, however, some experimental facts which should be carefully analysed in order to 
‚count for details of the theory. The most important among these seem to be the small 
syiations which have been observed in beta spectra by Langer's group (Johnson & al., 
amilton & al). The Fermi-Kurie plots for electron and positron spectra exhibited devia- 
ons from linearity which could not be explained by means of Fierz interference terms since, 
. present, there is strong experimental support (and even stronger theoretical reasons, too) 
r assuming these terms to vanish. Furthermore, introducing the Fierz interference terms 
ould imply opposite directions of the deviations for electron and positron spectra, whereas 
anger and his co-workers found experimentally that їп both cases there is an excess of 
w-energy ß-particles. In view of the somewhat strange consequences derived in this paper 
id following from the assumption of Langer's results, an experimental reinvestigation 
f the problem can prove to be necessary in order either to obtain support for our hypotheses 


r to show them to be wrong. 


'The following decays were investigated experimentally : 


BP 0. 325 (1+—> 0*), lg ft = 7.9, allowed G-T transition; 
En LP y 198p (1*— 07), Ig ft = 4.5, allowed G-T transition; 


* On leave of absence from the Institute for Nuclear Research, Warsaw, Poland. 
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3) 22Na > 22Ne (3+-> 24), lg ft = 7.4, allowed G-T transition; 
4) 90ү 907, (27 — 0%), lg ft = 8.0, first forbidden unique transition!. 
The analysis of experimental spectrum shapes resulted in the introduction of tl 


empirical multiplicative correction term: 


ix 


in order to obtain a linear Fermi-Kurie plot of the same type for В- and for ß*-spectra wit 
the values of b lying within the interval of 0.2 to 0.4 (in our system of units wherein =c- 
=т„=1). It was shown by the present author (Kuchowicz 1959a) that a spectrum shar 
of the desired type can be obtained under a simple assumption of semiphenomenologic: 
nature (addition of first order corrections in the gradient of lepton fields). The requiremen 
however, consists in fitting into a regular theory, without any ad hoc assumptions. 


II. Numerical estimations within the framework of the conventional theory 


It is well known that we have usually to deal in unique beta transitions with only on 
decisive energy-independent nuclear matrix element (Мот or Bj), which plays the rol 
of a multiplicative factor in constructing the Fermi-Kurie plot. With respect to the trans: 
lions mentioned, this is approximately true only in the case of the first forbidden decay c 
9 Y and in that of the allowed decay of !4In, whereas the two other allowed G-T decay 
exhibit too high values of lg ft, which is an argument in favour of the suppression of thi 
usually dominant matrix element. Thus, a series of additional correction terms usuall: 
considered irrelevant can be of comparable magnitude, even leading to a change of th 
resulting spectrum shape. In the conventional treatment of the (V — А) interaction, followin, 
additional terms for allowed G-T transitions occur: 


fen faxi [e [Go 


It was shown (Kuchowiez 1959b) that, in ordinary decays (in which lg ft does no 
exceed 4), these corrections are vanishingly small. This result can be adapted to the deca: 
of "In, whereas the two other decays (with matrix elements lesser by two orders of magni 
tude) escape this treatment and require a more thorough analysis. 

From the well known relation, valid for allowed transitions: 

СЕМеР + Cor Magi = 22102-0 a 


5 
t тос® 


we shall now estimate the Gamow-Teller matrix element (with value Cor = Са given by 
Goldhaber (1958)). This matrix element amounts to 7.10-3 for 32P and 12.5 - 10-3 for Na! 
We use the standard formula for the spectrum shape: 


еее» 


1 [n a paper published recently by Hamilton, Langer and Smith (Phys. Rev. 119, 772 (1960)) their result: 
have been extended upon the positron decay of 897г. ў 


TEE = Fy pg? Оса + ۱C4 | {| é 
| 
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x (25а е (Геза кв юе) 
«2 fa) (fa)! коз 2n (ССР CIC) O 
«el fs) (fuss) [5-25 }—зе|\дк 


and the approximate form of the nuclear matrix elements given by Morita (1959). According 
to the large ft-value, we assume that the matrix element, f © is approximately one percent 
of its value for ordinary favoured transitions. Thus, the ratios | Гу И o| etc. are multiplied 
by a factor of the order of 10?. Nevertheless, what we obtain is not simply a spectrum pro- 


j . b . . => > э,» > > 
portional to n 28 4 , since the terms in for?, [(a-7r)r and fa х r obviously imply energy 


dependence of a different type (terms proportional to E and E?) which failed to be observed 
hitherto. Only in one very special case do we obtain a spectrum of the desired empirical 
shape for ??P. This occurs if we take into account only the additional matrix element 
y ysr «0.0009 - 102/5; we obtain 50.29. Even under this rather artificial assumption of 
considering only one correction term no agreement with experiment can be achieved in 
the case of the positron decay of ??Na, since the coefficient b does not exceed 0.05 for the 
ratio |f/ygr|/|/o| taken from the reasonable interval of values (0+ 0.1). 


No explanation of the deviations in the beta spectra within the framework of convention- 
al corrections thus seems possible even in the case of the /-forbidden transitions. It seems ~ 
to be a mere accident that for one possible value of f ysr we obtained the empirical correction 
for 32P, since similar estimations fail to yield even approximately the correction factor for 
the positron decay of Na. There arises the problem: How to explain the deviations on 
a solid basis, not by assuming ad hoc unnatural cancellations of matrix elements, and without 
the necessity of introducing a different explanation for each experimental case. Besides, 
it is of interest to mention that there is one more phenomenon in beta decay that cannot 
be explained theoretically. Namely, experimental values of the longitudinal polarization 
of B-particles were found to be by about 10% lower than the theoretical value —v/c (Mikaelan 
and Spivak 1959). A reinvestigation of the problem by the same authors led them to the 
conclusion that the polarization had been systematically underestimated ; nevertheless, 
deviations from the theory, of the order of about 5% still remain. This problem seems to 
be closely related to the former one and, maybe, they might be solved simultaneously. 

In what follows, two tentative explanations of the deviations under consideration will 
be given, both of which are no more than the first step of a more detailed analysis to be 
carried out subsequently, and both requiring further support. One of these alternative 
explanations reduces deviations observed to possible G-nonconserving interactions, whereas 
the second will be seen to draw conclusions from the hypothesis of an intermediary chiral 


boson. 
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III. G-nonconserving -interactions 


In order to explain the deviations in the beta spectra, let us consider the most general 
: А m я X E cU 
transition matrix element for a nucleon decay їп Weinberg's notation (1958): 


M = [u,y,(1 + ys) uy] [us ГУ, + gyaj K" + thyk,) uy] + 
zr [u,iysy,(1 + ys) u,] ГАРАТА + га ys th 40; sk") u, | (3) 
The effect of strong interaction is described by means of the 6 form factors f 4, fy etc., which 
are functions of k?=k*k,, k, being the 4-momentum transfer. If time reversal holds and 
final state interactions are neglected, the form factors in (3) are real. 

Weinberg divided the primary weak interactions into two classes according to the 
G-transformation? properties of the strongly interacting currents. It was shown that, for 
nucleon processes, the terms with ^4 and hy can arise only from the second-class interactions 
which are absent in the Feynman-Gell-Mann theory. The presence of these interactions. 
would break the deep relation between the weak interactions and isotopie spin suggested 
by Feynman and Gell-Mann. We shall subsequently see that there is a test for the existence: 
of such interactions, 

Let us consider the beta decay of complex nuclei, assuming that the main contribution 
to the nuclear matrix elements comes from the structure of the nucleon rather than of the 
nucleus?, 

The relativistic form of the S-matrix element is given in Table 1. We apply the Foldy- 
Wouthuysen transformation, and the nucleon operators are taken in the nonrelativistic 
form as by Alaga &al. (1958). In Table I we have 4-component nucleon spinors, whereas: 


TABLE I 
Relativistic form of the S-matrix element for nucleon decay 
1 General form of the nucleon covariant for 
Lepton covariant - - 
pure Fermi transitions pure G-T transitions 
Я Д1) = [uz + у) u] lup(/v + gy P: a+ hy Е,В)и„]| lus(ey.P Вә) u,] 
Ца) = [uj (1 + уь) u] [us( — f y — ey Epa — hy. PB) uy] | [ug(— Јр gp Ey + igyP x Bo)u,l 
Lo) = [uz 0 (1 + ys) u,] [ug (ih, P x Ba) и, | lus Cf Ao gAPBy, — а Ева + 
A + ihgP x Ва) и] 
Му) = [uzys(1 + ys) u] осени lup(— Jays — gaEobys — ВАР: Ва) uy] 
КОЕР Е AT 
4 ОСЕ СЕЕ а = — —— 
Nonrelativistic form 2M 2M 
of the nucleon operators P 
ЮАР [02 
p=] ys = 
2M 
ЕИ 


The G-transformation is defined as the product of charge symmetry and charge conjugation. 
This means that we consider only the decay of a physical nucleon inside the nucleus, neglecting the 


possibility of exchange" or cooperative effects in complex nuclei. Since, for high Z, Coulomb effects destroy 
G-invariance, both first- and second-class interactions are retained. І 
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in Table II, where the nonrelativistic form of the S-matrix element is given, we have two- 
-component nucleon spinors and o is a Pauli matrix in nucleon space. 

In subsequent treatment, the dependence of the form factors on the four-momentum 
will be disregarded for a specific transition, these being treated as constants (which may, 
however, differ from one decay to another). A rough guess would yield the magnitude of 


TABLE II 


Nonrelativistic form of the S-matrix element for nucleon decay 


Lepton Nonrelativistic form of the nucleon covariant for 


covariant pure Fermi transitions pure G-T transitions 


2M 
: | EM 
A} | | om’ 84 TH ha) 220) P | u Í 
LG == E | ET +h diro E (Р -9) ША 
| (75) АГАЙ aa one 


The terms underlined are second order correction terms discarded in further cosiderations. 


the constants g and ^ at least about one order of magnitude lower than that of f. We shall 


= 


therefore disregard those terms in Table II wherein the g and A constants are multiplied 
by 2M (here M is the nucleon mass in units of electron mass). 


In order to obtain all desired formulae it is sufficient to adapt the results of Alaga’s 
paper after a suitable change of notation*. Alaga's terms with the nuclear matrix element До 
(for 1-forbidden transitions) are omitted in order to simplify the discussion. 


4 [n Alaga's formulae (1)... (9) we effect the following substitutions (with the new constants on the right 


запа side): 


a) for Fermi transitions: 


88 |fy + hy EP 


* 
| 1 
ue eee i hy Ey) Е 2 | 
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The problem of the relative signs for @—- and ß*-transitions may be easily solved i 
Weinberg’s Theorem I is taken into account. By this theorem, a change of the sign of л! 
and A, in the interference terms follows for f*-decay, and the final formulae for the spectrur 
shape are as follows: 


a) for Fermi transitions (0->0): 


Hy(E)dE = Fo(+ Z, E) pEg? 


ERE 2Re (fy Ay) Eo + Ви Е2 + 


Е 1 | 
+ 5 (и? 2Ке (fvhy) (1 + ED + 2Е, һә? | (а 2: — EI (4 


b) for Gamow-Teller transitions (AJ=]1): 


2 
M(E) dE = (+ 2, E) pEq*| L lua + 2Ке (faha) Eo + |ha|? E2 + 


Iur : : E : 2 SA ТЕ 
s 3 Е | fal? + 2Re (fahi) E a - 2E lhl | (ка 2Е — | + (E 


1 l مر رمم‎ Je 
#4: ane | yarar (ree) (m 28 #224 2) az 


The upper signs are for ß”-decay, the lower ones — for ß*-decay. The spectrum shapi 
for a AJ=0 (not 0>0) transition is given by the sum of expressions (4) and (5). By analog: 
to the formulae for the spectrum shapes, angular correlation formulae can be derived fror 
expressions (3) and (8) in Alaga's paper (1958). 

Langer's empirical shape correction factor is readily derived from (4) and (5). Let u 
assume for the time being, in order to simplify the respective expressions, that the constant 
fy ete. may be considered to be real numbers, and that the following estimation is correct 


1 
LAE z l 


1 1 
Ihal> Max = | /А|, M Ју], ivi (€ 


b) for Gamow-Teller transitions: 


84> lfa— Вай 


1 1 * 


2M 


Three other constants g,, g} and g, have to be put equal to zero. If the S, T and P interactions will come to 1 
considered once more, the formula presented by Alaga (1958) can be also adapted so as to take into ассош 
the additional terms for tensor interaction given by Weinberg (1958). We have applied here Alaga 


formulae with some sign corrections. It is the authors duty to thank Dr Tadié for sending him tl 
corrigendum to Alaga's paper. І 


1 1 1 * 
M 8384 2Re | — hgE)) [agros 2 | 
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It follows that a part of terms in (4) and (5) may be omitted and that the spectrum shape 
A : b 
is proportional to p Eq? ( + A , wherein 


a) for Fermi transitions: 


= 2хҥ(®ҥЁ + 1) 


Ду? Axp(&+ Ey) + х®Е (ЗЕ, + 48) 
* 
(p= 
Jv 
b) for Gamow-Teller transitions: 
bur E rr 2xer(xerky + 1) 
3+ 4xgT(É зе Eo) xerEo(Eo A£) 
хат = и 
Ла 


The upper (lower) signs are for 6-(87)-Чесау. 

The result thus obtained is important, though somewhat strange: There is some indica- 
tion that, in f-decay, we have an admixture of G-nonconserving interaction of the order 
of ten percent. For хст = — 0.255) we obtain 5— 0.42, 0.16 and 0.12 resp. for ??Na, 
32P and 140, which is consistent with the experimental values. This explanation of 

| Langer's results is still, however, no satisfactory argument in favour of the existence of 
` G-nonconserving interaction, and we must wait with final conclusions until more experi- 
mental data become available. Especially spectrum shapes for pure G-T and Fermi emitters 
should be carefully investigated in order to obtain further information on possible deviations. 


IV. The hypothesis of the intermediary singlet boson in beta decay 


In a series of papers by Tanikawa and co-workers (Tanikawa 1957, Tanikawa and 
Watanabe 1959, Oneda and Tanikawa 1959) a two-stage process is proposed for Fermi interac- 
tion, to the effect that the latter is transmitted by some Bose fields with zero spin having 
only renormalizable interactions with the Fermi particles. The Yukawa model of beta decay 
was modified so that the Bose field possessed a source consisting of a nucleon and a lepton. 
Similar modifications have been carried out with respect to other decays. The effective four- 
-fermion interaction in ß-decay is of the experimentally well established (V—A) form: 


Hg = вру 1+ ys) риф? (1 + ys) v, + Н.с. (7) 
'This interaction results by eliminating the B-field from the basic Hamiltonian: 
Hs = в 1 — ys) v£B + v, (1 — ys) v;B] + Н.с. (8) 


where B is a complex scalar neutral field with the same chirality as its factor in the Hamilto- 


nian Hç. y^ is the charge conjugate spinor Cy. 


5 The agreement of the theory with the experimental results for the ratio а Ла! Ма may be regarded 


as a post hoc justification of the estimation of eq. (6). 
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Some time ago various elegant methods and principles have been proposed for obtaining: 
the desired U.F.I. of (V—.4) form (Feynman and Gell-Mann 1958, Behrends 1958, Sakurai 
1958). None of these, however, could yield the experimental ratio of the axial vector and 
vector parts of the Universal Fermi Interaction. It was always — | instead of the experi- 
mental value — 1.25 for C4/Cy. One possible objection against the method proposed here: 
for introducing the (V— A) interaction is the same as the one to be raised against other: 
elegant principles from which this interaction can be deduced, viz.: it likewise leads to 
a C,/V ratio equal to — 1 instead of — 1.25. 

The intermediary boson hypothesis should not be regarded, however, as a mere artificial 
complication of the present situation. One argument in favour of this hypothesis is that! 
it explains in a straightforward manner the deviations in the beta spectra. In the energy, 
spectrum of the f — particles, in addition to the usual factor СР Ерд? (in standard B-decay 
theory notation), we obtain the factor 


Ф = (Mi — M2+2ME)2 (9) 


from the intermediate meson propagator. Here Mg. M and Е are the meson mass, nucleon: 
mass and electron energy, respectively. For a suitable value of the meson mass we can obtain | 
a spectrum shape of the desired type. There will be an excess of low energy f-particles 
of both signs and the deviations in the Fermi-Kurie plot may agree to a high degree of! 
accuracy with the deviations observed experimentally (Johnson 1958, Hamilton 1958) 
In order to illustrate this thesis we include a table of the values of 


а | 
3 p / 2 q (10) | 


which, on multiplication by = E,— E, should yield the Fermi-Kurie plot. The values are 
given for a f-transition in which the maximum f-particle energy E, is equal to 5 (in 
m,c* units). This value of E, is of the order appearing in Langer’s experiments (Johnson. 
1958, Hamilton 1958), where it amounted to e.g. 4.9 for "In and 4.34 for 32P, 

The numerical values in Table Ш are given with the accuracy to а common multiplica- 
tive factor chosen so as to yield unity for the spectrum end. 


TABLE IH 


Values of the quantity Э for a hypothetical f-decay with Еу-5 


Э from the formula with the 


А Э from the formula with the theoretical correction 

E empirical correction \- $i factor (9) for the boson mass Мр 

id b—0.2 50:3 1875 m, 1900 m, 1925 m, 1950 m, 
| 1.07 1.11 га 1.06 1.045 1.035 
2 1.03 1.04 1.08 1.045 1.035 1.025 
3 1.02 1.02 1.05 1.03 1.02 1.015 
4 1.005 10015 1.025 1.02 1.01 1.01 
5 | | 1 l 1 1 f l 
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It is reality concluded from Table III that the possible existence of an intermediary 
boson of mass ~1900 m, is able to account satisfactorily for the deviations in beta spectra 
observed. 


V. Concluding remarks 


In proposing the foregoing alternative explanations for the deviations in beta spectra, 
the author did not intend to exhaust all the theoretical possibilities, but aimed at confining 
attention to those only which seemed the simplest in the present stage, when the scarcity 
of experimental data provides no possibility of distinguishing between them. | 

The author wishes to express his gratitude to prof. I. А. Smorodinsky for his enlightening 
discussions of the subject. | 
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GAMMA-VIBRATIONAL LEVELS IN Er 
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( Received November 10, 1960) 


Directional correlation measurements in 188Ег for the gamma-gamma cascade 779.8— 
80.59 keV were performed. The values of spin /—3 and the quantum number K— 2 for the energy 
level 860.2 keV could be determined. The deformation of the 166Er nucleus in the ground state 
and in the gamma-vibrational state, as well as the rotational-vibrational interaction of energy 
levels in this nucleus, are discussed. 


1. Introduction 


In the course of the last year some measurements were performed [1—6], on the basis 
of which level schemes for 16Er in the decay of 1867 were proposed. One group of authors 
[1, 3, 6] based their proposals on measurements of energies and partly also intensities of 
conversion electrons, another group [2, 4, 5] — on gamma-gamma coincidence experiments. 
All these experiments are consistent with the presence in 166Er of a ground state rotational 
band with excited levels at about 80 keV (0,24), 265 keV (0,4+) and 545 keV (0,67). However, 
in the case of the next group of levels, the agreement of different authors in interpreting 
them is not so good. Boskma and De Waard [2] proposed а gamma-vibrational level (2,2) 
at 862 keV with a first rotational state (2,3) at 957 keV. The other authors [1, 3, 4, 5, 6] 
suggest the gamma-vibrational level at 788 keV and two rotational states at 862 keV (2,37) 
and 957 keV (2,44). It is very important to know whether there exists in 199Ег a gamma- 
-vibrational level and if so, to know the exact positions of the levels. These informations 
would allow us to say something about the deformation of the nucleus in the excited states 
as well as about the rotational-vibrational interaction. The experimental branching ratios 
of transitions from the levels in question to those belonging to the ground state band fit 
better the proposal of a y-vibrational band with the band head at 788 keV and two rotational 
states above it [4]. The accuracy of the measurements is however too low to confirm unambi- | 
guously this interpretation. The accuracy of the determination of the conversion coefficients 
is also too low to give the exact multipole order of transitions to the ground state band. 


* Joint Institute for Nuclear Research, Dubna, SSSR. 
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One can only say that the experimental results do not exclude the possibility of this transi. 
tions to be of the E2 type. 

The investigations reported in this paper were undertaken in order to obtain the spin 
and А values of at least one of the levels 788, 862 or 957 keV by means of the gamma-gamma 


directional. correlation measurements. 


2. Experimental 


For the angular correlation measurements in 198Ег ytterbium sources were used which 
were obtained by irradiation of Ta targets with 660 MeV protons in the synchrocyclotron 
of the Joint Institute for Nuclear Research in Dubna (USSR). The chromatograpbically 
seperated Yb fraction was then purified twice in the chromatographic column. The ytterbium 
lactate was evaporated on an aluminium foil. The diameter of the source was less than 
1.5 mm. The measurements in Cracow started two days after the irradiation. 

Two scintillation counters with RCA 6810 4 photomultipliers and NaJ(T1) crystals 
13" x 1" were used. The resolution time of the fast-slow coincidence circuit 2r = 3.75 x: 

1078 sec. 

The correction factors for the finite angular resolution of both counters were determined 
experimentally, before the correlation measurements started, applying the method of Lawson: 
and Frauenfelder [7]. For this purpose gamma rays of the Ni cascade were used. The width! 
of the collimated gamma ray beam was measured with a NaJ(T1) crystal 2 mm x2 mm x: 

«7 mm. The whole width of the beam at a 410 mm distance from the end of the collimator 
was 7,2 mm. The angular efficiencies g(a) and eH(x) for both scintilation counters were: 
measured at a distance of 76 mm between the source and the fronial surfaces of the crystals, 
the two photopeaks of Ni being selected with a single channel analyzer. The integrals: 
1a, ГР (cos a)e(a)|sin &|d were evaluated numerically. So obtained values of Iz 1g =0.9701 
and of /,/1, —0.900 differ from the theoretical ones, for the same distance and the y-ray energy: 
1250 keV [8] only by 0.3% and 1.794 respectively. These deviations are smaller than thei 
experimental errors, thus for the solid angle correction in the correlation experiments in 
166Кг we used the theoretical coefficients given by Stanford and Rivers [8] 

The gamma ray spectrum in the decay Yb (58h)—>166Tm (7.7h)—>188Er is very compli- 


cated and the half-time is rather not very long. Therefore, we were forced to chose such: 


. 


а у-у cascade for the spin determination of levels in question for which a good selection 
of the transitions is possible and also the differences in the anisotropy values for different: 
spins of the initial state are as large as possible. The table I gives the anisotropy values for: 
у-у cascades in which the first y-ray corresponds to the transitions starting from the: 
786.7, 860.2 or 957.0 keV levels, and the second y-ray corresponds to transitions between 
the levels of the ground state rotational band in 166Ey nucleus (fig. 1). Assuming, that the 
upper level band is connected with y-vibrations (К = 2+) and that the selection rules: 
for the quantum number K are exactly fulfilled, we can accept that both transitions in the 
cascades are of the E2 type. In table I only the values of spins proposed by Boskma and De: 
Waard [2], as well as the others [1, 3, 4, 5] are taken into account. One can see from this: 


table, that the cascades best suited for the determination of the spins of the initial levels 
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y correlation measurements are the following ones: 876.3— 80.59 keV, 594.9—184.5 keV 
nd 779.8—80.59 keV. Unfortunately, a sufficiently good selection of the 876.3 and 594.9 keV 


ransitions is by scintillation methods impossible, thus only the cascade 779.8— 80.59 keV 
suitable for measurements. 


15453 (0 2+) 


4617 0 0+ 
O 
a | | a 
N N 
1246.1 


54539 0 6 + 


166 Er Е (KeV) 
Fig. 1 


The correlation measurements for the 779.8—80.59 keV cascade were performed at 
| distance of 73 mm between the source and the frontal faces of the crystals. Care was taken 
when discriminating the y-rays corresponding to the transition 779.8 keV to select this 
is good as possible from y-rays corresponding to transitions 158.9 keV (tripple cascade 
/58.9—876.3—80.59 keV) and 876.3 keV (cascade 876.3 —80.59 keV) (fig. 1)-„Аз thick 


TABLE I 
Theoretical values of anisotropy A for cascades I;(E2)1(E2)];. 


Energy of 
transitions 


Н 692.1— 184.5 876.3— 80.59 594.9— 184.5 119.8— 80.59 706.0— 80.59 
in cascade 


(keV) 


Spin se- 
quence A—4—2| 3—4—2 | 4—2—0| 3—2—0|3—4—2 | 2 4—2 | 3—2—0 2—2—0/2—2—0| — 


I; 


Anisotropy | 743 | —930 | +16.6 | —33.3 | —9.30 | +38.3 | —333 | +7.69 | +7.69| — 
A(%) 
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absorber Pb-Cu-Cd protected the crystal in order to cut off the low energetic y-rays. Or 
the erystal of the second counter, selecting the 80.59 keV transition, another Pb-Cu-Cc 
absorber was placed, the thickness of which was selected experimentally to give a strong 
reduction of the X-ray peak and only a small weakening of the 80 keV photopeak. А narrow 
channel width of the channel analyzer was set on the top of the 80 keV photopeak. The 
source was centered uniil the single counting rates in both detectors for different positiona 
of the moveable counter were constant to less than 0.5%. Every 10 minutes the positiom 
of one counter was changed from 90° to 180° and back, and coincidences as well as single 
counting rates in both channels were recorded. 


The following anisotropy of the gamma-gamma cascade 779.8—80.59 keV was obtained 
А = — (20.9 + 2.8)% 


The exact correction of this result for the coincidences of overlapping cascades is 
impossible because of the complicated y-spectrum in the 166Yb decay. Nevertheless, the 
approximate calculations based on coincidence and single counting rates, setting the channels 
of the analyzers on the peaks and near by the peaks, show that the result being obtained: 
may be considered to be lowered by less than 1%. 


3. Discussion 


After the correction for the finite angular resolution in our experiments, the theoretical 


779.8 keV 80.59 кеу... € 
> 2+ > 0+ in !#Er are A = — 29.80% für 


I; = 3 and A = + 4.13% for J; = 2. Our result points rather the 3—2—0 sequence, although 
the experimental value is too low. However, the intermediate state 80.59 keV has a measured. 
half-life T, = 1.7 x 1059 sec [9], and the 166Er nucl 


anisotroples for the cascade 1, 


eus is strongly deformed. The quadrupole 
interaction in the intermediate state will thus have a strong smearing out effect on the correla- 
tion pattern. This effect accounts partly for the deviations between the theoretical and 
experimental anisotropies. 

Marklund et al. [10] have measured у-у directional correlations in !99Er from the 
"Ho (27h) decay, with the same intermediate state 80.59 keV 


. They used aqueous solu- 
tions of holmium trichloride as sources. 


These authors give the values for the atten- 
uation factors G,=0.63 +0.03 and G,=0.69 + 0.03. Assuming the same attenuation in 
our case the anisotropy in the cascade 779.8—80.59 keV for 1,= 3 is expected to be A= 
=—19.49%. This value is in agreement with our result within the experimental errors. 

Up to now we have assumed that the 779.8 keV transition is a pure £2 radiation. Let 
us now discuss the possible admixture of М1 in this transition. Without taking into account 


the attenuation of the correlation caused by the finite angular resolution and by the quadru- 


pole interaction in the intermediate state, we obtain for the anisotropy the value 4 = —23.31% 
for the initial spin I; —2 and a 20% admixture of M1 radiation and 6<0 in the 779.8 keV 
transition. If the 860.2 keV level is connected with y-vibrations (K=2+), this high admixture 
of M1 radiation is rather improbable. If the spin of the 860.2 keV level 1]; —3, then a 1% 


admixture of M1 gives A= 


A=— 41.19% for 6>0 and A=—22.38% for д<0. In the first 


case the attenuation has to be very strong, in the second сазе — neglectable to fit the experi- 
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ental result. Larger admixtures of М1 are giving larger deviations from the theoretical 
alue A= —33.35% for pure E2 transitions. Thus we can accept that the 779.8 keV transi- 
on is a pure £2 one. A possible admixture of M1 radiation is less than 1%. 


. E 779.8 keV 80.59 keV 
The anisotropy measurements of the cascade [,———2* rn 


0+ point that the 
60.2 keV level has a spin 7, —3 and K —2. That means, it is the first excited level of the 
otational band connected with the y-vibrational level 786.7 keV 22: 

Fig. 1 represents the incomplete level scheme of 196Ег. The new energy values of levels 
jere calculated on the basis of mean values of transition energies given by Jacob et al. [3] 
nd Grigoriev et al. [6]. This level scheme contains three excited rotational levels in the 
round state band and three members of the gamma-vibrational band. 

In the 16Ho (27h) decay a beta-vibrational level (0,0) in 188Ег at 1460.2 keV was 
eported [10]. The level at 1461.7 keV excited in the decay of 1887 [3, 5, 6], may be identical 
ith the level (0,07). It is also possible that the level 1545.3keV [4, 5, 6] is another rotational 
tate (0,2+) connected with /-vibrations. 

On the basis of the proposed level scheme, and using the rotational energy formula 


аа AT Tele KOKA DONE B [TOC END) UG]? (1) 


here E, is the energy of the band head, 4 =h?/2J with J — moment of inertia of the nucleus, 
— the rotational-vibrational interaction coefficient, we are able to calculate the constants 
and В. for the ground state band and the beta and gamma vibrational bands. 

The hydrodynamical approximation gives the rotational-vibrational coefficient for 


ЛЛ ЛЬ 4 
ва - (27) FE ma] ©) 


‚here ho, and ho, are the excitation energies of beta and gamma vibrations respectively. 

The calculated values of A and B for the ground state band and for beta and gamma 
ibrational bands, as well as the factor b>=B,/B,,, for the ground state band are shown 
n table II. From this table it is seen that 4 —?/2/ is somewhat larger for the ground state 


die 


he ground state band: 


TABLE II 


Moments of inertia and rotational-vibrational interactions in '®*Er nucleus. 


n DIO. 
= — (keV В. (eV) De 
A 27 ( е ) ехр же 
sround state ground state ground state 
е 5 В х am 
band band an 
13.50 | 1229 | 1.90* 12.4 | 9.5 2.39 0.14 


* Calculated with B=0 


and than for the y-vibrational band. If the simple formula J = 3B,P*, connecting the mo- 
rent of inertia J with the deformation В, is approximately valid, then one can say that the 
SEr nucleus is still more deformed in the excited state than in the ground one. 
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The ratio Б= Bal Bexp for the ground state band in 188Ег is somewhat larger than th 


value b=2.12 proposed by Sheline [11] to be a constant for all nuclei. | 
The authors are indebted to Mr. Chojnacki and Mr. Yu. Norseev for the work don: 
with the preparation of sources. Thanks are also due to Miss Z. Konieczna, Miss M. Miesowic 


and Mr. G. Zapalski for their valuable assistance in the course of measurements. 
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PHASENSPRUNG IM BRENNPUNKT 


Von A. RUBINOWICZ 
Polnische Akademie der Wissenschaften 
( Eingegangen am 28. November 1960) 


Im Jahre 1938 hat der Verfasser auf Grund der Kirchhoffschen Theorie der Beugung und 
unter teilweiser Benutzung geometrischer Methoden bewiesen, dass der Phasensprung im Brenn- 
punkt eine geometrisch-optische Erscheinung ist. Er tritt nämlich bereits in der nach den Ge- 
setzen der geometrischen Optik konstruierten einfallenden Lichtwelle auf. Es wird nun für diese 


Tatsache ein mit Hilfe rein analytischer Methoden durchgeführter Beweis angegeben und es 
werden weitere Plausibilitätsbetrachtungen über die Entstehung des Phasensprunges im Brenn- 
punkte mitgeteilt. 


In einer vor längerer Zeit veröflentlichten Abhandlung hat der Verfasser (Rubino- 
wicz 1938) gezeigt, daß der Phasensprung im Brennpunkt ein geometrisch-optisches Phüno- 
men ist. Er ist nämlich bereits in der nach den Gesetzen der geometrischen Optik zu defi- 
nierenden einfallenden Lichtwelle vorhanden. Der Beweis dieser Behauptung wurde unter 
Zugrundelegung der Kirchhoffschen Beugungstheorie durchgeführt. In. Übereinstimmung 
mit den Youngschen Anschauungen über die Entstehung der Beugungserscheinungen 
wurde nämlich die durch das entsprechende Kirchhoffsche Integral gegebene Wellenbewe- 
gung in eine im Sinne der geometrischen Optik definierte einfallende Lichtwelle und eine 
von dem beugenden Rande ausgehende Beugungswelle aufgespalten. Dies wurde in Anleh- 
nung an eine früher vom Verfasser (Rubinowicz 1917) vorgenommene Zerlegung der Wel- 
lenbewegung durchgeführt, die im Falle einer divergierenden Kugelwelle durch das Kirch- 
hoffsche Integral dargestellt wird. Dieses Verfahren erfordert aber kompliziertere geo- 
metrische Betrachtungen!, die nicht jedermanns Sache sind, so daß es gerechtfertigt 
erscheinen mag, einen mit Hilfe der im allgemeinen gebräuchlicheren analytischen Methoden 
durchgeführten Beweis der damaligen Ergebnisse mitzuteilen. Wir tun dies nun nach dem 
Muster eines später vom Verfasser (Rubinowicz 1953) für eine divergierende Kugelwelle 
angegebenen analytischen. Verfahrens zur Zerlegung der durch das Kirchhoffsche Integral 
gegebenen Wellenbewegung in eine einfallende und eine Beugungswelle. Im Anschluß 


1 Die geometrischen Betrachtungen, die in meiner Arbeit über die Beugung einer divergierenden Kugel- 
welle (Rubinowiez 1917) sowie in der Arbeit über den Phasensprung (Rubinowicz 1938) enthalten sind, können 
mit Hilfe der Vektorrechnung etwas einfacher dargestellt werden. Vgl. die in Vorbereitung befindliche englische 
Ausgabe der „Beugungswelle“ (Rubinowiez 1957). 
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an diese im Sinne der Youngschen Vorstellungen vorgenommene Aufspaltung der Weller 
bewegung sollen auch neue Plausibilitätsbetrachtungen über die Entstehung des Phaser 


sprunges im Brennpunkte angegeben werden. 


§ 1. Das Vektorfeld des Kirchhoffschen Ansatzes und sein Vektorpotential 


Sei f eine Fläche, die eine durch den beugenden Rand B begrenzte beugende Öffnun 
in einem Schirme 5 verschließt (vgl. Abb. 1). Durch / hindurch möge eine konvergierend: 
Lichtwelle mit dem Brennpunkt F einfallen. Die auf diese Weise entstehende Wellenbewe 
gung wird in der Kirchhoffschen Beugungstheorie (vgl. etwa Rubinowiez 1957) auf de 
Seite des beugenden Schirmes, auf der sich der Brennpunkt F befindet, durch das Integre 


1 | eir e e— iko е- iko e eikr | 
6 =Z 4 1.1 
ще Ar. | г on 0 о дп|г !| af ( 
da 


- 


AES : 3% a | 
beschrieben. Dabei wird hier der Zeitfaktor exp (—iwt) vorausgesetzt. о = ЕО bzw. r= PC 
bedeutet die Entfernung eines Punktes О auf der die Beugungsöflnung verschlieBender 
Fläche f vom Brennpunkte F bzw. vom Beobachtungspunkte P. 


ol 


Abb. 1. Zweidimensionale schematische Darstellung der Beugung einer konvergierenden Kugelwelle. S ist dei 
Schirm, B der beugende Rand, J eine die Beugungsóffnung verschliessende Fläche, F der Brennpunkt, P eir 
Beobachtungspunkt, О ein Punkt auf der Fläche J; К, bzw. K, sind die Schattengrenzen des vorderen bzw. 
rückwärtigen Lichtkegels. Durch + bzw. — wird die positive bzw. negative Seite der Schattengrenze bezeichnet. 


Um die Wellenbewegung (1.1) im Sinne der Youngschen Beugungstheorie umzuformen. 
bemerken wir, daß der Integrand in (1.1) durch die Normalkomponente des Vektorfeldes 


x 1 e tke eiko Л » 
V(Q) = |: grad —— = ДЕ] à Í (1.2) 


0 
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jauf der Flüche f gegeben wird, wenn wir es als Funktion des Punktes Q auffassen, der nun 
und in den zunächst folgenden Überlegungen irgend einen Raumpunkt bedeuten soll. 
Da e*"/r und ео als Funktionen von Q aufgefaßt Lösungen der Schwingungsgleichung 
Zu --k*u —0 sind, so ist div V=0. Es muß sich daher das Vektorfeld V durch ein Vektor- 
potential W darstellen lassen, so daß 


V = го W. (1.3) 


ist. Gelingt es nun das Vektorpotential W anzugeben, so ist die Wellenbewegung WE’), CE) 
in der Gestalt 
u(P) — ү. rot, Wdf 
if 

darstellbar. Die Umformung des Flächenintegrals (1.1) für (Р) in ein die Beugungswelle 
darstellendes Randintegral iiber den beugenden Rand B wird dann unmittelbar durch den 
Stokesschen Integralsatz gegeben. 

Um nun das Vektorpotential W zu berechnen, bemerken wir, daB mit Rücksicht auf 
grad o = о/о und grad r = r/r das Vektorfeld V (1.2) nun die Gestalt 


З ee Л AN 
ааа (е | Aed e (14) 
erhält. > 


Da die Vektoren 0 und F ihren Beginn in den Punkten F und P haben, so müssen die 


Vektoren V des Vektorfeldes (1.4) stets in einer durch diese beiden Punkte gehenden 
Ebene Е liegen. Das Vektorpotential W darf demnach keine in der Ebene EFegenden Kompo- 
nenten besitzen. Würde dies nämlich der Fall sein, so müßte das aus (1.3) berechnete Vektor- 
feld V nach der Definition der Rotation auch in den zu den Ebenen E senkrechten Rich- 
tungen nicht verschwindende Komponenten aufweisen. Wir dürfen daher für das Vek- 


torpotential W den Ansatz 
W(Q) =r x e &(Q) (1.5) 


machen, wo g(Q) eine skalare Funktion des Punktes Q ist. 
Mit Rücksicht auf die Beziehung rot(r x o) — 2(r — o) folgt dann aus (1.3) 


V = g rot (r x 2) + (grad g) x (r x 9) = 2(r — 0) g+ rl grad g) — lr grad g). (1.6) 
Beachten wir, daß 
1 3 
2g+ r grad g = = r grad (r?g) 


ist, und eine analoge Beziehung mit o und o statt rund г besteht, so können wir (1.6) auch 
st, б 0 


in der Gestalt 


= 


pam nest (o grad (о?) a (r grad (r?g)) 


darstellen. Setzen wir diesen Ausdruck für V gleich dem Ausdrucke (1.4) und multiplizier en 
wir die erhaltene Beziehung vektoriell einmal mit r und das andere Mal mit о, so ergeben 
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sich schlieBlich für die Funktion g die beiden Differentialgleichungen 


ik (r— o) 1 | 

ik(r—o) : |, | 

X 1 Е : E 1) = o grad (020). (1.7Ъ` 
«кл 


Um zunächst eine Lösung der ersten Differentialgleichung d. h. (1.7a) zu finden spalter 
wir von der in (1.7a) auftretenden Funktion rg den Exponentialfaktor exp ik(r—o) ab. 
setzen also 


pag = giMr-Og (1.8) 


Diese Maßnahme wird durch die Tatsache nahegelegt, daß ein Exponentialfaktor bei allen 
Differentiationen erhalten bleibt und daher (1.7a) durch den Ansatz (1.8) von dem Expo- 
nentialfaktor exp ik(r—o) befreit wird. 


Mit Rücksicht auf 
1 => 
r grad (r — o) = r1 [ro — ro] 


sowie den Ansatz (1.8) erhalten wir für die rechte Seite der Diflerentialgleichung (1.7a) den 


Ausdruck 


T grad r?g = eik(r-o x [ro — rà] f + ейс-9 r grad f. (1.9). 


- 


Setzt man (1.9) in (1.7a) ein, so tritt auf beiden Seiten der so erhaltenen Beziehung 


LT YE 1 tk 
In l> +a] = 0 TOILET grad f (1.10) 
4ло\о о 0 

Фе Konstante Ё nur einmal und zwar nur in der ersten Potenz auf. Nimmt man nun an, 


daß in der Funktion f die Konstante k nicht mehr enthalten ist, so erhält man die beiden 
Beziehungen 


* 


DW 1 
f- Am o ro — Fe’ (1.11a) 
+ 17 
r grad f = ee (1.11b) 


Unter den gemachten Voraussetzungen wird somit die Funktion J durch den Ausdruck 
(1.11a) gegeben, muß jedoch auch noch die Beziehung (1.11b) erfüllen. Dies ist in der Tat 
der Fall, wie man mit Hilfe der aus grad ro=r+9 folgenden Relation 


Sor 0 ПА Зали cepa 
r grad Е [ro — ro] = = [ro — ro] 


bestätigt. 
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Aus (1.8) und (1.11а) ergibt sich schließlich, daß die Funktion 


1 eik(r—e) 1 
En An ro re—ro 


(1.12) 


ine Lösung der Differentialgleichung (1.7a) ist. 

Dieser Ausdruck für g muß auch noch die Differentialgleichung (1.7b) befriedigen. 

ies ist in der Tat der Fall, weil ja (1.7a) bei Vertauschung von r, r und o, o sowie Vorzeichen- 
änderung von г in (1.7b) übergeht, während sich dabei die durch (1.12) gegebene Funktion g 
icht andert. 
- Zu den obigen Ausführungen sei noch bemerkt, daß man den gleichen Ausdruck 
1.11а) für die Funktion f erhält, wenn man voraussetzt, daß f zwar von k abhängt, jedoch 
in eine nach fallenden Potenzen von Ё fortschreitende Potenzreihe entwickelbar ist. Setzt 
an nümlich? 


f= AT (1.13) 


in (1.10) ein, so erhält man für die Funktionen f, die Beziehungen 


E TET 
| Genel Suma. gr (1.14a) 
| г - TOT 
e DU TO. iue r grad fo, (1.14b) 
re rol fex = — T grad. fa (n z 1). (1.140) 
0 


Aus (1.14a) folgt, daß fo die durch (1.11a) gegebene Gestalt hat. Da infolgedessen f, die 
Relation (1.11b) erfüllt, so ergibt sich aus (1.14b) daß f, = 0 ist. Setzt man f, = 0 in die 
Rekursionsformel (1.14c) für n —1 ein, so folgt daß auch f,=0 ist. Benützt man die gleiche 
Relation (1.14c) der Reihe nach für n=2, 3, 4, ... so kann man schließen, daB auch alle 
übrigen f, verschwinden und daher die Funktion f (1.13) durch (1.11а) gegeben wird. 


$ 2. Die Beugungswelle 


Aus (1.3) und (1.5) folgt, daB die durch den Kirchhoffschen Ansatz (1.1) gegebene 
Wellenbewegung auch in der Gestalt 


u(P) — i rot, (r x og) df (2.1) 
f 


larstellbar ist, wobei g durch (1.12) gegeben wird. 


? Der gleiche Ansatz (1.13) kann auch im Falle der Beugung einer divergierenden Kugelwelle (Rubino- 
vicz 1953, 1957) zur Auflösung der dort auftretenden zu (1.11a) und (1.11b) analogen Gleichungen verwendet 


verden. 
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Bei der Anwendung des Stokesschen Satzes auf (2.1) sind nun die beiden Fälle zv 
unterscheiden: 

(a) In keinem Punkte der die beugende Öffnung verschließenden Fläche f wird da: 
durch (1.5) und (1.12) definierte Vektorpotential W unendlich. 

(b) Es ist ein solcher singulärer Punkt Q, auf der Fläche f vorhanden. 

Um festzustellen in welchen Fällen solche singuläre Punkte Q, auftreten, bemerke: 
wir: Damit W unendlich wird muß der Nenner ro — ro in g (1.12) verschwinden, d.h. e 
muß ro = ro sein. Das kann nur eintreten, wenn die beiden Vektoren r und 0 die gleich! 
Richtung und den gleichen Sinn haben. Es verschwindet jedoch dann auch der Vektor r x ¢ 
Um das Verhalten von W in einem solchen Falle zu untersuchen bezeichnen wir mit В de: 


Winkel zwischen den beiden Vektoren * und o. so dab 


Ir x o| sin f 
2 ro | — cos В 


= сіб б (2.2 


wird. Falls jedoch die beiden Vektoren und о die gleiche Richtung und den gleichen Sin: 
haben, so wird f—0 und daher der Ausdruck (2.2) unendlich. 

Der Fall paralleler Vektoren und 0 kann nur eintreten (vgl. Abb. 1), wenn der Beobach 
tungspunkt P eine solche Lage hat, daß die durch ihn und durch den Brennpunkt Е hin 
durchgehende Gerade die Fläche f in einem Punkte trifft, der wie oben mit Qo bezeichne 
werden möge. Beobachtungspunkte P, denen solche singuläre Punkte Q, auf der Fläche‘ 
entsprechen, befinden sich somit auf Geraden, die durch den Brennpunkt F hindurchgehe 
und die Flüche f schneiden. Die Gesamtheit aller dieser Geraden liegt jedoch innerhal 
des Doppelkegels, dessen Scheitel mit dem Brennpunkt F zusammenfällt, und desse 
Mantel aus allen durch F und den beugenden Rand B hindurchgehenden Geraden besteh! 
Diesen Doppelkegel, in dem alle Beobachtungspunkte liegen, denen singulire Punkte Ç 
auf der die Beugungsóffnung verschließenden Fläche f entsprechen, wollen wir als de 
Lichtkegel bezeichnen. Seine Mantelfläche soll die Schattengrenze heißen. Wie sich nämlie 
in 83 herausstellen wird, tritt in dem Lichtkegel auch eine Wellenbewegung auf, die wi 
als die einfallende Lichtwelle bezeichnen kónnen. Der Halbkegel, dessen Basis durch di 
Fläche / gegeben wird, soll im folgenden der vordere Lichtkegel, der andere Halbkegel de 
rückwärtige Lichtkegel heißen. 

Wir beschäftigen uns zunächst mit dem Falle (a), wo der Beobachtungspunkt P au: 
serhalb des Lichtkegels liegt. Da dann der Integrand in (2.1) keinen singulären Punkt Оу, а 
der Fläche f besitzt, liefert die Anwendung des Stokesschen Satzes auf (2.1) nur ein Integr 
über den beugenden Rand B des Schirmes S, welches mit Rücksicht auf (1.12) gegebe 
wird durch 


EUIS CE 
pra mem то то 
В 


(2. 


Dieser Ausdruck stellt eine von dem beugenden Rande В ausgehende Wellenbewegur 
dar und ist daher im Sinne der Youngschen Vorstellungen über die Entstehung der Be 
gungserscheinungen als die Beugungswelle anzusprechen. Aus der oben gemachten Vorausse 
zung (a) folgt, daß der Beobachtungspunkt P außerhalb des Lichtkegels liegt, also in eine 
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Raumgebiet, das als der Schattenraum bezeichnet werden kann. Bei der Annüherung des 
Beobachtungspunktes Р an einen Punkt P, der Schattengrenze wird der Integrand in der 
Beugungswelle ив (2.3) unendlich. Und zwar findet dieses Unendlichwerden für jenes Randele- 
ment ds, statt, das auf der Verbindungsgeraden des Brennpunktes F mit dem Punkte 
Po liegt. Für den ds, auf dem beugenden Rande entsprechenden Punkt werden nämlich 
die beiden Vektoren * und 0 einander parallel. Dies hat aber mit Rücksicht auf (2.2) ein 
Unendlichwerden des Integranden in (2.3) zur Folge. Diese Tatsache bewirkt daf) die 
Beugungswelle ug (2.3) in der Schattengrenze nicht stetig ist, nümlich hier einen Sprung 
besitzt, der, wie im nachfolgenden Paragraphen klar werden wird, notwendig ist um den 
entgegengesetzt gleichen Sprung der einfallenden Welle zu kompensieren. 

Bezüglich der weiteren Diskussion des Ausdruckes (2.3) für die Beugungswelle ив 
verweisen wir auf die früheren diesbezüglichen Publikationen des Verfassers (Rubinowicz 
1938, 1957). Hier findet der Leser auch die entsprechenden Näherungsformeln für die 
Beugungswelle (2.3). 


$3. Einfallende Lichtwelle 


Liegt der Beobachtungspunkt P innerhalb des vorderen oder rückwärtigen Lichtkegels, 
so schneidet die durch Р und den Brennpunkt F hindurchgehende Gerade die Flüche f, so 
daß hier ein singulärer Punkt (y des Vektorpotentials W auftritt. Wir haben es dann mit 
dem Falle (b) des $2 zu tun. Bei der Anwendung des Stokesschen Satzes auf (2.1) müssen 
wir daher das Kurvenintegral nicht nur über den beugenden Rand B sondern auch über 
eine den singulären Punkt 0, auf der Fläche f umgebende Kurve А erstrecken. Während 
jedoch der Umlaufssinn des Kurvenintegrals über den beugenden Rand B mit der Nor- 
malenrichtung auf der Flüche f durch die Rechtsschraubenregel gekoppelt ist, wird der 
Zusammenhang zwischen dem Umlaufssinn auf der Kurve K und der Normalenrichtung 
auf dem aus der Fläche f durch die Kurve K herausgeschnittenen Flächenstück durch die 
Linksschraubenregel gegeben. 

Das Kurvenintegral über den beugenden Rand B ergibt ebenso wie im Falle (a) die 
Beugungswelle ug (2.3). Das nun zusätzlich auftretende Integral über K stellt hingegen die 
einfallende Lichtwelle und zwar sowohl im vorderen als auch im rückwärtigen Lichtkegel 
dar. Um dieses Integral über K zu berechnen, wollen wir auf А statt der Bogenlünge s eine 
neue Integrationsvariable а einführen. Zu diesem Zweck betrachten wir die Schar aller 
Halbebenen, die durch die durch die beiden Punkte P und Е hindurchgehende Gerade 
begrenzt werden und bezeichnen mit x den Winkel zwischen einer beliebigen Halbebene £ 
und einer fixen Halbebene Ед dieser Schar. 

Um nun den Zusammenhang zwischen da und dem Randelement ds auf dem beugenden 
Rande anzugeben, benennen wir mit R die Entfernung der beiden Punkte P und F und 
mit d den Abstand eines Randelementes ds der Kurve K von der durch die beiden Punkte P 
und F hindurchgehenden Geraden. Da4 Rd den Flächeninhalt des Dreieckes F, P, ds 


bedeutet, so stellt 


r тхо 3.1) 
Ар ( 
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einen Einheitsvektor dar, der senkrecht auf der Ebene des Dreieckes F. Р, ds steht. Па 
Nds = ds, die Projektion des Randelementes ds auf die Richtung des Einheitsvektors M 
bedeutet, so ist dæ = ds,/d, so daß mit Rücksicht auf (3.1) 


(r x o) ds = + Ка?ах (3.2) 


ist. Wir haben die rechte Seite von (3.2) mit dem + Zeichen versehen, da wir die Integra- 
tion über & zwischen den Grenzen 0 und 2л erstrecken wollen. Es gilt dann in (3.2) das 
positive bzw. das negative Vorzeichen, je nachdem der Winkel y zwischen dem Vektor 
r x 0 und dem Bogenelement ds auf der Kurve K ein Spitzer oder ein stumpfer ist. Falls 
nur die Kurve K eine solche Gestalt hat, daß der Winkel x beim Durchlaufen, der Rand- 
kurve sich stets nur monoton ündert, kann man leicht aus der Abb. 1 entnehmen, dab у 
einen spitzen bzw. einen stumpfen Winkel darstellt, je nachdem sich der Beobachtungspunkt. 
P im vorderen bzw. im rückwürtigen Lichtkegel befindet. Demnach ist in (3.2) das positive 
bzw. das negative Vorzeichen zu verwenden, je nachdem der Beobachtungspunkt P im 
vorderen bzw. im rückwürtigen Lichtkegel enthalten ist. 

Den Nenner in g (1.12) kann man mit Hilfe der nachstehenden, durch Quadratbildung: 
aus R = o — r folgenden Beziehung 


2[ro — rg] = R? — (o — r)? (3.3) 
ausdrücken. 
Um noch die in (3.2) auftretende Entfernung d durch R, г, o darzustellen, bemerken 
wir, daß aus der Heronschen Formel für den Flächeninhalt des Dreieckes F, P, ds sich die 
Beziehung 


Rid? = 1 [(r+ of — К] [R? — (r — oy] (3.4) 


ergibt, da ja seine Seiten durch R, r, o gegeben werden. 
Bei der Anwendung des Stokesschen Satzes auf (2.1) erhält man daher für das über 
die Kurve А erstreckte Integral mit Rücksicht auf (3.2), (3.3) und (3.4) den Ausdruck 
2z 


1 eik(r—e) ^ А 
u 7i [(r + o)? — А] da. (3.5) 


c. 
0 
Dabei gilt mit Rücksicht auf (3.2) im vorderen Lichtkegel das positive und im rückwärtigen 
das negative Vorzeichen. 
Bezeichnen wir mit ry bzw. оу die Absolutbeträge der nach dem Punkte Qo der Fläche у 
7 > zur е 4 "1 = 1 1 
von P bzw. F aus gerichteten Vektoren, so erhalten wir beim Zusammenschrumpfen der 
Kurve K auf den Punkt О, im Grenzfalle für (r+ o)? — R? = (r+ 0)? — (To — 00)? den Wert 
Ато. Der Abstand R des Beobachtungspunktes P vom Brennpunkt F wird nämlich im 
vorderen Lichtkegel durch 
К = og — re (3.68) 
und im rückwärtigen Lichtkegel durch 


R= To — 00 , (3.65 


gegeben. 
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Es ergibt sich daher im Grenzfalle des Zusammenschrumpfens für den Beitrag (3.5) 
des Kurvenintegrals über K zu der durch den Kirchhoffschen Ansatz (1.1) dargestellten 
)Wellenbewegung im vorderen bzw. im rückwärtigen Lichtkegel mit Rücksicht auf (3.62) 
bzw. (3.6b) für die einfallende Lichtwelle der Ausdruck 


e-ikR 
ив, = <. (3.72) 
i bzw. 
eikR 
ug, = — =. (3.7b) 


jim vorderen Lichtkegel wird somit die einfallende Lichtwelle (3.7a) durch eine auf den 
| Brennpunkt F hin konvergierende Kugelwelle gegeben, wie man dies gemäß der geometri- 
schen Optik erwartet. Im rückwärtigen Lichtkegel verdient (3.7b) eigentlich nicht die 
Bezeichnung „einfallende“ Lichtwelle, da sie ja doch eine von dem Brennpunkt F aus 
] divergierende Kugelwelle darstellt. Wenn wir sie so benennen, so geschieht es nur um für 
die beiden durch die geometrische Optik gegebenen Wellenbewegungen ug, und ug, eine 
gemeinsame Bezeichnung zur Verfügung zu haben. 

Die einfallende Welle ug, (3.7a) im vorderen Lichtkegel unterscheidet sich von der 
einfallenden Welle ug. (3.7b) im rückwärtigen Lichtkegel durch das Vorzeichen. Es ist somit 
ein Phasensprung z bereits in der einfallenden Lichtwelle vorhanden. Damit ist unsere 
Behauptung bewiesen, daß der Phasensprung im Brennpunkt eine geometrisch-optische 
Erscheinung ist. 

Allerdings kann man dagegen geltend machen, daß der durch den Kirchhoffschen 
Ansatz (1.1) beschriebene Wellenvorgang im ganzen Raume also auch im Brennpunkte F re- 
gulär ist. Das bedeutet, daß die in der einfallenden Welle ug, + ug, im Brennpunkte F 
auftretende Singularität durch die Beugungswelle ug kompensiert werden muß, also im 
Brennpunkte F selbst gar nicht zur Geltung kommt. Und zwar gilt diese Tatsache ohne 
Rücksicht auf die Größe der Wellenlänge, bleibt also auch im Grenzfalle 4>0 erhalten. 

Diese letztere Tatsache wird dadurch bedingt, daD im Brennpunkte F die Beziehung 
r=o gilt (vgl. Abb. 1, in der jedoch anzunehmen ist, daß der Punkt Qo auf dem beugenden 
Rande В liegt) und daher gemäß (2.3) von allen Punkten des beugenden Randes die 
Beugungswelle ohne Rücksicht auf die Größe von A mit der gleichen Phase eintrifft. In 
einem Beobachtungspunkte P, der ganz in der Nähe des Brennpunktes F liegt, ändern sich 
demnach die Phasen der von den einzelnen Punktes des beugenden Randes eintreffenden 
Beugungswellen nur wenig, so daD die gesamte Beugungswelle ив hier noch gréBere, d. h. mit 
der einfallenden Lichtwelle vergleichbare Werte erreichen kann. Sobald aber A gegen Null 
strebt, «o werden die Phasen der von den einzelnen Punkten des beugenden Randes kommen- 
den Beugungswellen sich immer rascher andern und daher immer mehr auslóschen (vgl. 
z. B. Rubinowicz 1957, S. 166 ff.). Das bedeutet jedoch, daß der Bereich um den Brennpunkt, 
‘in dem die Beugungswelle mit der einfallenden Welle vergleichbare Werte annimmt und 
daher die Wellenbewegung von der durch die einfallende Lichtwelle gegebenen betrüchtlich 

abweicht, zugleich mit A immer kleiner wird. Dementsprechend wüchst der Bereich in dem 
| die geometrisch-optische einfallende Lichtwelle als eine angenäherte Darstellung der gesamten 
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Wellenbewegung angesehen werden kann, in dem also der Phasensprung sich geltend macht... 
In diesem Sinne ist somit die Behauptung zu verstehen, daß der Phasensprung eine вео-. 


metrisch-optische Erscheinung ist. 


$ 4. Plausibilitätsbetrachtungen über die Entstehung des Phasensprunges 


Nachdem vom Standpunkte der Kirchhoffschen Beugungstheorie der mathematisch 
strenge. Nachweis für die geometrisch-optische Natur des Phasensprunges erbracht ist,. 
drängt sich die Frage nach seinem intuitiven Erfassen von selbst auf. Ich habe zwar ini 
meiner früheren diesbezüglichen Arbeit (Rubinowicz 1938) mit Hilfe von Plausibilitits- - 
betrachtungen eine Begründung für unsere These gegeben; sie hat jedoch den Nachteil,. 
daß sie nicht voll der Tatsache gerecht wird, daß die Beugungswelle ug im Brennpunkte F 
eine Singularität aufweist. 

Um Plausibilitätsbetrachtungen durchzuführen, die von diesem Mangel frei sind,. 
gehen wir von der Bemerkung aus, daß durch den Brennpunkt die Schattengrenzen der 
von allen Punkten des beugenden Randes stammenden Anteile der Beugungswelle ug, 
hindurchgehen. Nähern wir uns einem Punkte der Schattengrenze, der nicht in unmittelbarer 
Nühe des Brennpunktes liegt, so bewirkt dies, wie wir es in § 2 auseinandergesetzt haben, 
für einen bestimmten Anteil der Beugungswelle ein Unendlichwerden des Integranden. | 
Nähern wir uns nun dem Brennpunkt, so nähern wir uns gleichzeitig den Schattengrenzen 
der von allen Punkten des beugenden Randes stammenden Anteile, wodurch also das 
Auftreten von Singularitäten in allen diesen Anteilen bewirkt wird. Auf diese Weise wird 
es verständlich, daß die Beugungswelle im Brennpunkte so stark unendlich wird, daß sie 
das Unendlichwerden der einfallenden Welle ир, + ug, kompensiert. 

Wir wollen noch zeigen, wie man auch verschiedene andere Details des Verhaltens der 
Beugungswelle in der Nachbarschaft des Brennpunktes verstehen kann. Sie muß nämlich 
zunächst die im vorderen und rückwärtigen Lichtkegel auftretenden und hier verschiedene 
Vorzeichen aufweisenden Singularitäten (3.7a) und (3.7b) der einfallenden Welle beheben. 
Außerhalb dieser beiden Lichtkegel muß sie jedoch auch in der Nachbarschaft des Brenn- 
punktes endlich bleiben. Den Schlüssel zum Verständnis dieses Tatbestandes liefert das 
Verhalten der Beugungswelle wenn wir uns einem Punkte ihrer Schattengrenze von der 
einen oder der anderen Seite nähern. Die Grenzwerte, die wir in diesen beiden Fällen er- 
halten, unterscheiden sich nämlich in erster Näherung nur durch das Vorzeichen, haben 
aber sonst gleiche Absolutbeträge und gleiche Phasen. Der Nachweis dafür kann in dem 
betrachteten Falle ganz ebenso wie im Falle der divergierenden Kugelwelle (Rubinowicz 
1957, S. 194) erbracht werden. Das Auftreten. dieses Sprunges der Beugungswelle wird 
durch die in § 2 besprochene Singularität des Integranden von ug (2.3) für dasjenige Randele- 
ment ds, bewirkt, das auf der Verbindungsgeraden des Brennpunktes F und des Punktes IS 
auf der Schattengrenze liegt, dem der Beobachtungspunkt P sich nähert. Die Verschiedenheit 
des Vorzeichens rührt von dem Anteil x 0) ds des Integranden in (2.3) her, der zum 
Randelement ds, gehört. In der Schattengrenze selbst, wo r und о parallel werden, verschwin- 
det nämlich (r x o)ds und ändert beim Durchgang durch die Schattengrenze sein Vorzeichen. 
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Der Bequemlichkeit des sprachlichen Ausdruckes wegen, bezeichnen wir nun diejenige 
eite der Schattengrenze eines Randelementes ds als die positive Seite, die vondem vorderen 
ichtkegel aus zu erreichen ist. Die andere Seite heiBe die negative. Wie aus der Abb. 1 zu 
ntnehmen ist, náhern wir uns nun der positiven bzw. der negativen Seite der Schatten- 
renze, wenn wir an sie von dem vorderen bzw. rückwürtigen Lichtkegel aus herantreten. 

Bei der Annäherung an den Brennpunkt F vom vorderen bzw. rückwärtigen Lichtkegel 
us, wird daher das Verhalten der Beugungswelle ug durch die Grenzwerte auf der positiven 
zw. auf der negativen Seite der Schattengrenze bestimmt (vgl. Abb. 1). Dies erklürt die 
erschiedenheit der Vorzeichen der Beugungswelle in den beiden Lichtkegeln, die notwendig 
t, um die Singularitüten der einfallenden Wellen (3.7a) und (3.7b) im Brennpunkte Ё zu 
eheben. Die Tatsache, daß der Grenzwert der Beugungswelle endlich ist, wenn wir uns 
em Brennpunkte Е aus dem Außenraume der beiden Lichtkegel nähern, erklärt sich 
araus, daD wir uns dann gleichzeitig der negativen Seite des vorderen und der positiven 
eite des rückwärtigen Lichtkegels nähern (vgl. Abb. 1). Bei der Annäherung an den Brenn- 
unkt F von dem Aufenraume der beiden Lichtkegel her, haben somit die Singularitüten 
erschiedene Vorzeichen und ergeben daher einen endlichen Grenzwert. 

Dabei ist noch zu beachten, dab die Beugungswelle ug (2.3) aus Wellen besteht, die ihre 
usgangspunkte in dem beugenden Rande haben. Sie kann daher nur singulüre Wellen 
rzeugen, die diese Fortpflanzungsrichtung haben, die also in dem vorderen Lichtkegel auf 
en Brennpunkt F hin, im rückwürtigen aber von ihm weg gerichtet sind, genau so wie dies 
ei den einfallenden Wellen (3.7a) und (3.7b) der Fall ist, deren Singularitüten durch die 
eugungswelle beseitigt werden. 
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